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Abstract 

We present Lyapunov stability and asymptotic stability theorems for steady state solutions of 
general state-dependent delay differential equations (DDEs) using Lyapunov-Razumikhin meth¬ 
ods. Our results apply to DDEs with multiple discrete state-dependent delays, which may be 
nonautonomous for the Lyapunov stability result, but autonomous (or periodically forced) for 
the asymptotic stability result. Our main technique is to replace the DDE by a nonautonomous 
ordinary differential equation (ODE) where the delayed terms become source terms in the ODE. 
The asymptotic stability result and its proof are entirely new, and based on a contradiction ar¬ 
gument together with the Arzela-Ascoli theorem. This approach alleviates the need to construct 
auxiliary functions to ensure the asymptotic contraction, which is a feature of all other Lyapunov- 
Razumikhin asymptotic stability results of which we are aware. 

We apply our results to a state-dependent model equation which includes Hayes equation as a 
special case, to directly establish asymptotic stability in parts of the stability domain along with 
lower bounds on the size of the basin of attraction. 
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1. Introduction 

We consider the following general delay differential equation (DDE) 
discrete state-dependent delays, 

u{t) = f(t, u(t), u(t - T\(t, u(t))), . .., u(jt — r,v(t, u(t)))), 
u(t) = 

and prove Lyapunov stability and asymptotic stability results using Lyapunov-Razumikhin tech¬ 
niques. We apply our results to the model state-dependent DDE 

ii(t) = fiu{t) + cru(t - a — cu(t)), t ^ 0 , 

u(t) = (p(t), 0, (L2) 


in d dimensions with N 


t > t 0 , 
t ^ to. 


( 1 . 1 ) 
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with a > 0, as an example of (1.1), to directly show asymptotic stability in parts of the stability 
domain, and derive bounds on the basin of attraction. 

Differential equations with state-dependent delays arise in many applications including milling 
[21], control theory [41], haematopoiesis [6] and economics [29]. There is a well-established 
theory for retarded functional differential equations (RFDEs) as infinite-dimensional dynamical 
systems on function spaces [15,7, 10], which encompasses problems with constant or prescribed 
delay, but very little of this theory is directly applicable to state-dependent delay problems. Ex¬ 
tending the theory to state-dependent DDEs, including equations of the form (1.1) is the subject 
of ongoing study. See [16] for further examples and a review of recent progress. 

The model state-dependent DDE (1.2) includes the constant delay DDE often known as Hayes 
equation [17] as a special case when c = 0. Hayes equation is a standard model problem used to 
illustrate stability theory for constant delay DDEs in most texts on the subject including [15, 39, 
20], as well as being a standard numerical analysis test problem [2]. Hayes equation is also used 
to illustrate Lyapunov-Razumikhin stability results in [1, 15, 37], 

The state-dependent DDE (1.2) was introduced by Mallet-Paret and Nussbaum, and is a natural 
generalisation of Hayes equation to a state-dependent DDE with a single delay which is linearly 
state-dependent. But whereas Hayes equation is linear, the DDE (1.2) is nonlinear and can admit 
bounded periodic solutions. Mallet-Paret and Nussbaum investigate the existence and form of 
the slowly oscillating periodic solutions of a singularly perturbed version of (1.2) in detail in [35] 
and use it as an illustrative example for more general problems in [32, 33, 34]. This DDE ia also 
studied in [3, 18, 19, 24,31], 

Following [15], to define an RFDE let be the ^-dimensional linear vector space over the 
real numbers equipped with the Euclidean inner product ■ and Euclidean norm | • |. Let r ^ 0 and 
C = C([-r, 0], W 1 ) be the Banach space of continuous functions mapping [—r, 0] to R' / with the 
supremum norm denoted || • ||. If u e C([?o - r, tf], then for every t e [to, tf] define u t e C by 

u,{6) = u(t + 9), 9e[-r, 0]. (1.3) 

Then for F:lxC-> M. d , with the dot denoting a right-derivative, an RFDE is defined by 

u(t) = F(t, u t ), u to = <p e C. (1-4) 

A solution of the RFDE (1.4) is a function u 6 C([to - r, t/],R rf ) which satisfies (1.4) for t e 
[f 0 - r, tf). If F is Lipschitz onKxC then existence and uniqueness of solutions is assured [15], 
and the RFDE (1.4) defines a dynamical system with the function space C as its phase space. 

If a constant function <p e C yields Fit. ip) - 0 for all t then ip is a steady state of the RFDE. To 
study the stability of steady states of RFDEs, the method of Lyapunov functions for ODEs was 
first extended to Lyapunov functionals V : RxC-» 1 for RFDEs by Krasovskii [25], Lyapunov 
theorems for stability of RFDEs require the time derivative of the functional along a solution of 
(1.4) to be nonpositive or strictly negative, similar to the theorems for ODEs using Lyapunov 
functions [15, 25]. However, finding functionals k:RxC->R with this property for RFDEs is 
much harder than in the ODE case where the Lyapunov functions have the form V : R x W 1 —> R. 

Razumikhin [38] developed the theory on how one might go from using the more difficult 
Lyapunov functionals for RFDEs back to Lyapunov functions again. His fundamental idea is 
that it is only necessary to require a constraint on the derivative of V whenever the solution is 
about to exit a ball centered at the steady state. Following this approach, Barnea [1] presents 
a Lyapunov stability theorem for RFDEs and also considers Hayes equation. A comprehensive 
discussion of Lyapunov functionals and functions for general RFDEs is presented by Hale and 

2 


Verduyn Lunel in chapter 5 of [15]. Other works with Razumikhin-type results include [14, 22, 
23, 25, 26, 27, 37, 43]. Of these [27, 37, 43] include results tailored for time-dependent delays. 

State-dependent DDEs of the form (1.1) can easily be written as RFDEs. For example, with 

Fit, ip) = pcp(O) + (Tipi-a - cip(Q)), (1.5) 

the model problem (1.2) is an RFDE of the form (1.4). However, two problems arise. Firstly, we 
need an a priori bound r on the delays in (1.1) to consider it as an RFDE. This can be overcome 
for specific problems such as (1.2) for which bounds can be derived under certain parameter 
conditions. The second difficulty is more fundamental; it is well known that state-dependency of 
the delays results in F not being Fipschitz onlxC [16]. This is easily shown directly for (1.5). 

The Fyapunov-Razumikhin results of Barnea [1] only establish Fyapunov stability and assume 
that F is Fipschitz and so are not applicable to state-dependent DDEs of the form (1.1) when they 
are written as RFDEs. Other authors, such as Hale and Verduyn Funel [15], make the weaker 
assumptions on F, but use auxiliary functions to establish Fyapunov stability and uniform asymp¬ 
totic stability. The construction of these functions is nontrivial in all but the simplest examples, 
and we have not seen such functions constructed for a state-dependent problem. Rather than try 
to circumvent these problems for RFDEs, in Section 2 we will develop new proofs of Fyapunov 
stability and asymptotic stability for the state-dependent DDE (1.1) with N discrete delays. 

In Section 2 we present our main Fyapunov-Razumikhin stability results. In Assumption 2.1 
we state the assumptions that we make on the nonautonomous DDE (1.1) with N (state- 
dependent) delays, the main ones being that / is locally Fipschitz with respect to its arguments in 
M. d and the delays are locally bounded near to the steady state. Then in Theorem 2.5 we provide 
sufficient conditions for Fyapunov stability of a steady state of the DDE. The main idea behind 
the proof is the conversion of the DDE into an auxiliary ODE problem where the delayed terms 
are regarded as source terms. In Theorem 2.7 we establish asymptotic stability of the steady state 
when the DDE (1.1) is autonomous. For simplicity of exposition we present the proof for the 
case of a single delay, but the result remains true for multiple delays or periodically forced prob¬ 
lems. This result is significantly different to previous Fyapunov-Razumikhin asymptotic stability 
results which require auxiliary functions to establish uniform asymptotic stability. In contrast. 
Theorem 2.7 does not require the construction of any auxiliary functions, and is proved using the 
auxiliary ODE by a contradiction argument, which shows there cannot exist a solution which is 
not asymptotic to the steady state. 

Theorems 2.5 and 2.7 establish Fyapunov stability and asymptotic stability when the solutions 
of the auxiliary ODE have certain properties, but to determine those properties exactly would 
require the solutions of the DDE. So, in Section 2, we also show how to define a family of 
ODE problems that are subject to constraints defined by bounds on the DDE solution and its 
derivatives, which can be determined without solving the DDE. Lemma 2.3 establishes bounds on 
the growth of solutions to ( 1.1 ) which is used to ensure solutions remain bounded for sufficiently 
long (k times the largest delay for some integer k) to acquire k bounded derivatives. Stability is 
then established from the solution properties of this family of constrained ODE problems. 

In Section 3 we review the stability region of the model equation (1.2) which is known [12] to 
be the same for the state-dependent (c 4 0) and constant delay cases (c = 0). We also consider 
the properties of the auxiliary ODE we define for this problem, and define sets and functions that 
are required in the following sections to apply our Lyapunov-Razumikhin results. 

In Section 4 we apply Theorem 2.7 to provide a proof of asymptotic stability of the steady state 
of model equation (1.2) in subsets of the known stability region, together with lower bounds on 
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the size of the basin of attraction. This result is given as Theorem 4.7. Since the delayed inputs 
to the auxiliary ODE have k — 1 bounded derivatives, for the k — 2 and k = 3 results we establish 
suitable bounds on the first and second derivatives of the ODE source terms, while the k - 1 
result, does not require any differentiability of these terms. 

The expressions for the stability regions derived in Sections 4 all involve a term that needs to be 
maximized over a closed interval (see the definitions of the P(6, c, k ) functions in Definition 3.5). 
For k = 1 this maximum is readily evaluated, while for k - 2 an expression for the maximum 
is established in Theorem 4.2 (whose proof is given in Appendix A). Plots and measurements of 
the derived asymptotic stability regions in (p, cr ) parameter space are given in Section 5, where it 
is seen that these regions grow with the integer k, but do not appear to fill out the entire stability 
region in the case p + 0. In Section 5, we also briefly review previous work on the c - 0 constant 
delay case of (1.2) with p — 0 and p 4- 0 and point out an error in the results of Barnea [1]. 

In Section 6 we present two examples of solutions which are not asymptotic to the steady state 
of the model equation (1.2) when p ^ 0 > cr, and which hence give upper bounds on the radius 
of the largest ball contained in the basin of attraction of the steady state. We compare these with 
the lower bounds on the basin of attraction given by Theorem 4.7 for k - 1,2 and 3. In Section 7 
we present brief conclusions, and compare and contrast our approach with linearization. 

2. Lyapunov-Razumikhin techniques for state-dependent DDEs 

Here we state and prove our main theorems to establish the Lyapunov stability and asymptotic 
stability of steady state solutions to state-dependent DDEs of the form (1.1). 

We will consider continuous initial functions ip e C, where for the state-dependent DDE (1.1) 
we let C — C([-r(d), 0], R rf ), with r(d) defined by (2.1) below. By a solution of (1.1) we mean a 
function u e C'([to, tf), M . d ) which satisfies (1.1) for t e [to, tf) for tf > to and u(to) = if (to)- We 
denote by B( 0, d) the closed ball centred at zero with radius d in R d . We will assume that (1.1) 
has a steady state at u = 0 and make the following assumptions throughout. 

Assumption 2.1. Let d, N and keZ, d^l, N^l, k^l and to e R. 

1. / : R x Rdv+tjd an{ j T . . |- fo , R for i = 1,..., N are continuous functions 

of their variables. 

2. f(t, 0,0,..., 0) = 0 for all t ^ to- 

3. T;(f, 0) R 0 for all t ^ to and i = 1,..., N and the constant r max = max, = i ^ sup f>fo r,(t, 0) 
satisfies T max e (0, oo). 

4. There exist positive constants Lq, L\,..., and 5o such that 

\f(t,u,Vi,...,v N ) - f(t,u,v i,...,Vat)| < Lo\u -u\ + Li|vi - vi| + • •• + L N \v N - v N \ 

for all 1 ^ to and u,v i,..., vn, m, vi, ..., v# e B( 0, do) c R 0 '. Let L - Lq + L\ + • • • + Ln. 

5. The delay terms T,(t, u ) are nonnegative and Lipschitz continuous in u on [to, °°) X B((), do) 
with Lipschitz constants L Tj . 

6. f(t, u, Vi. \’n) is at least maxj/: — 2, 0} times differentiable in its u and v variables, and 

t ft, u) is at least max{A' —1,0} times differentiable in u. 

Definition 2.2. For any d e (0, do] define 

r(d) = sup r i(t, u). 

i=l,...,N 
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( 2 . 1 ) 


Assumption 2.1, items 3 and 5 imply that 0 < r max ^ r(S) f T max + S max, L Ti < oo with 
lim 0 .. ( o rid) = r max > 0. In contrast, the bound on the delay terms used in RFDE theory by Hale 
and Verduyn Lunel [ 15] is a global bound on r,(f, u) over all t and it. 

The Lipschitz continuity conditions in Assumption 2.1 , items 4-5 are based on those of Driver 
[8], From the results in [8], these conditions ensure local existence and uniqueness of a solution 
to (1.1) close to the steady state solution for all sufficiently small Lipschitz continuous initial 
functions ip, and existence of a solution for all sufficiently small continuous initial functions ip. 

We will prove Lyapunov stability and asymptotic stability by contradiction, by showing that 
there is no solution of (1.1) that does not have the required stability property. Thus we do not 
need to make assumptions specifically to ensure existence or uniqueness of solutions. Rather, we 
make sufficient assumptions so solutions that do exist have the properties we require. 

Since with state-dependent delays it can be hard ensure a priori that the delays r,(t, u(t)) are 
strictly positive along solutions, and problems with vanishing delays can be interesting, we do 
not assume that delays T;(f, u) are strictly positive in Assumption 2.1. This makes our results 
more widely applicable, allowing us to deal with state-dependent DDEs without having to first 
impose or prove that the deiays do not vanish aiong solutions. The possibility of vanishing delays 
will on the other hand complicate some of our proofs, starting with Lemma 2.3 below. Similar 
results to Lemma 2.3 are well known for time-varying and nonvanishing delays and proved using 
a Gronwall lemma and the method of steps; however the method of steps is not applicable with 
vanishing delays, and we will need a more technical proof. 

The following lemma provides bounds on the growth of solutions to ( 1. 1). For initial functions 
ip € C we will use Part II of Lemma 2.3 to ensure that solutions u(t) remain bounded for a 
sufficiently long time interval so that u(t ) acquires the differentiability that we require. Part I of 
Lemma 2.3 will be essential in the proof of asymptotic stability in Theorem 2.7. 

Lemma 2.3. Suppose that Assumption 2.1 is satisfied. 

I. Let T > to and e > 0. Suppose that u(t) and u(t) are solutions to (1.1) with initial functions 
ip and ip e C respectively, where both |<^(f)l ^ So and \ip{t)\ ^ So for t ^ to, and \u(t)\ ^ So 
and |m(0| ^ So for all t e [to, T], If also \ ip(t) — ip(t)\ ^ efor t ^ to then 

\u(t) - u(t)\ ^ ee L(f - to) Vr e [f 0 , T], 

II. Let T > to and 5 e (0, ho]- Let \<p(t)\ ^ Se L( -'°~ T) for t ^ to- Then any solution of (1.1) 
defined for t 6 [to, T\ satisfies 

\u(t)\^Se U ' , - T) Vte[f 0 ,r], 

Proof. This proof is similar to the derivation given in Halanay [13] for RFDEs of the form (1.4), 
and which Halanay credits to Krasovskii. We prove part I first. For L > L suppose there exists 
t, e [to, T) such that Part f holds with L replaced by L only up to t - t t . That is 

| u(t) - S(t)| ^ 

for all t e [to, f»], but that for every t > t, there exists te (f„, t] such that \u(t)-u(t)\ > The 

simplest situation would be where for some h > 0 we have |«(t)-M(t)| > seL'~'^ for t e (f„, t»+/j), 
but the existence of such an interval is neither guaranteed by the conditions of the lemma, nor 
needed in this proof. First we note that 

I u(t) - u(t)\j t \u(t) - u(t) | = jj t \u(t) - u(t) | 2 = [u(t) - u(t)] ■ [u(t ) - u(t)] < \u(t) - u(t) | | u(t) - ii(t)\. 
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Since |n(f„) - u(t*)\ = ee L(f * to) > 0, dividing by the common factor leads to 


I HU) ~ u(U) | ^ f t \u(t) - u(t)\\ t=h . 


( 2 . 2 ) 


Since the solutions to (1.1) have right derivatives for t (p to, we can also compute 



Since r* - r,(r«, < t* it follows that \u(t„ - t,(7„, w(f»))) - n(f» - r,(L, w(f«)))| ^ ee L< -‘- But, 

using Assumption 2.1, item 4 we see that 


|«(f«) - S( f*)| = | /(f», w(f,), n(t* - Ti(f«, M(f«))),..., w(f» - Tjv(r», w(f„)))) 


- /(f», w(L), u(t„ - ri(f„ u(t*))),..., S(f» - Tjv(r„ S(r«))))|, 
< Lo|w(f«) - fi(f,)| + Li|t<(r* - Ti(t«, n(f.))) - u(t, - n(f., fi(f.)))| 


... + Lat|m(?* - TAr(f., w(t«))) - f<(f* - Tat(?*, n(f»)))|, 
^ (Lo + L x + ... + L N )ee L{, '-' 0) = Lse U, ’-' o) < Lee I( '*“ ?o) . 


Thus \ii(t*) - u{t,)\ < j t \u(t) - «(0l| f _, , but this contradicts (2.2), hence there does not exist such 

a f„, and so \u(t) - u(t)\ ^ ee^ (Mo) for all t e [to, T], Since this holds for all L> L, part I follows. 

We prove part II by contradiction. If part II is false then for some initial function ip such that 
\<p(f)\ ^ (5e L(f °~ r) for t ^ to there exists T* e (to, T) and a solution u(t) such that \u(T*)\ > 6e L{T ~ T) 
but \u(t)\ ^ 5o for all t 6 [fo, T*]. Letting ip(t) = 0 for all t < 0 from Assumption 2.1, item 2 we 
also have the solution u(t) = 0 for all t ^ to- Now applying part I on the interval [to, 7’* ] with 
these initial conditions and s = 6e L(u> , ] we find that 


\u(T*)\ = \u(T*) - u(T*)\ < Se L(, °- T) e L{r '- ,o) = 6e L(T ’- T) , 
which contradicts \u(T*)\ > de L,T ‘ 71 and so part II of the lemma holds. 


□ 


The growth bounds on the solution in Lemma 2.3 are valid even in the case of vanishing delays. 
If Tj(t, u) were bounded away from zero, a Gronwall argument could be used to obtain a tighter 
bound, but we allow for the possibility of vanishing delays in the state-dependent case. 

In this section we will develop a constructive technique for showing Lyapunov stability and 
asymptotic stability based on Lyapunov-Razumikhin ideas. To establish Lyapunov stability in 
Theorem 2.5 we will show that solutions remain in a closed ball B( 0, 6) of radius 6 about the 
steady state. We are thus implicitly using the Lyapunov function V(u) - [m| 2 /2, though it will not 
appear directly in our results. Asymptotic stability is established in Theorem 2.7, by showing 
that all solutions that remain in B( 0, 5) must converge to the steady state. The application of these 
results to the model problem (1.2) is demonstrated in Sections 3-4. 

We will prove Lyapunov stability (and later asymptotic stability) by contradiction. Suppose 
the DDE (1.1) has a solution which escapes the ball B( 0, 6) at time t, then u(t) = x with \x\ = 6. 
Now letting v(0) = u(t + 6) and 


77,(0) = u(t + 6 — Tj(t + 9, u(t + 6))), 6 e [~r(<5), 0], i = 1,... ,N 

for this solution u(t), we can rewrite the DDE ( 1. 1) as a nonautonomous ODE 


m = f(t + 6, vie), 77i(0),..., ij N ie)) 
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where 


v(-r t(t, x)) - uit - Ti(t, x)) - 77 /( 0 ), i — 1,..., N. 

Hence for each i such that r,(f, x) > 0 the escaping solution of the DDE (1.1) corresponds a 
solution of an ODE boundary value problem (BVP) with v(-r,(f, x)) - 77 ,( 0 ) and v(0) = x. To 
establish stability it is sufficient to show that the ODE BVP does not have any such solutions. We 
will use Lemma 2.3 part II to ensure that \u(t)\ xj d for t SC to + kr{6) so that the forcing functions 
77 ,( 0 ) acquire the regularity we require. Hence we define the set of forcing functions 77 , which 
could correspond to an escaping trajectory as follows. 

Definition 2.4. Suppose that Assumption 2.1 is satisfied for (1.1) and k ^ 1. Let 8 £ (0, bo], 
\x\ = 8 and t to + kr(8). Define the set. 


E (k) (8, x, t) 


: 77 / e C k 1 ([-t-(c5),0],.B(0,c 5)), 
and conditions 1 and 2 are satisfied. 


(2-3) 


1 . x ■ f(t, x, 771 ( 0 ),..., ? 7 /v( 0 )) ^ 0 , 

2. For some initial function ip £ C, equation (1.1) has solution u(t) such that 77 ,( 0 ) = u(t + 6 - 
Tj(t, u(t + 0)) for 0 € [-r(d), 0] for each ie {1,..., N }. 

Condition (1) in the definition is equivalent to ^|w(r)l ^ 0, a necessary condition for the 
solution to escape the ball B( 0, 8) at time t. In the following theorem we prove the first of our 
main results, that if certain conditions hold for all functions in the sets E^)(8, x, t ) then the steady 
state of (2.9) is Lyapunov stable. The condition (2.5) implies that the ODE BVP discussed above 
has no solutions, while ( 2 . 6 ) allows solutions with ^|v( 0 )|e=o ^ 0 . 

The sets E^(8, x, t) however cannot be determined without solving (1.1), so it is not practical 
to actually solve for them. Instead, the conditions of the theorems can be shown to hold for larger 
sets that contain E^){8, x, t). We prove the theorems first and then consider such larger sets. 

Theorem 2.5. Suppose that Assumption 2.1 is satisfied for (1.1). Ford e (0,do], x e \x\ — 8, 
define E^(8, x, t) as in Definition 2.4. Consider the family of auxiliary ODE problems, 

I v(0) = fit + 0, v(0), 771 (0),..., ij N (0)), 0 e [-T,-(r,x),0], 

\ v(-Ti(r,*)) = 77,(0) £R d , ^ ; 

for i = 1,... ,1V and t to + kr{8). We denote the solution of (2.4) by v(x, 771 ,..., if 

we want to emphasize the dependence on x and r/i, or just v(0) otherwise. Suppose there exists 
<5 1 6 ( 0 , do] such that for all 8 £ ( 0 , hi), and for every x such that |jc| = 8 and all t to + kr(8), 
for all ( 771 ,..., 77 /v) £ E( k) {8, x, t) the solution of (2.4) for some I £ {1,... ,N] satisfies rfit, x) > 0 
and either 

|v(x, 771 ,..., ?7jv)(0) ■ x < 8, (2.5) 

or 

|v(x, 771 ,..., 77 /vXO) • x = 8, and v(x, 771 ,..., 77 /vXO) ■ x ^ 0 , ( 2 . 6 ) 

then the zero solution to (1.1) is Lyapunov stable. Moreover if 8 £ (0,di) and |^(s)| < 8e^ Lkrt ^ 
for s £ [to — 7 "(d), fo] then the solution of ( 1 . 1 ) satisfies \u(f)\ ^ 8 for all t to. 

Proof. Let the hypothesis of the theorem hold and let 8 £ (0, d)). For A e (0, I) suppose |^>(s)[ ^ 
A8e~ Lkr (® for v £ [fo-K<5X foL By Lemma 2.3 part II, the solution of (1.1) satisfies u(f) £ B( 0, A 8 ) 
for all t SC to + kr{8). We will prove that nit) € B{ 0, d) for all t f to by contradiction. 
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Figure 1: Illustration of the proof of Theorem 2.5 in one dimension. 


Assume that a solution escapes the ball B( 0, 6 ) for the first time at some time t* > to + kr(6), 
then there exists 6 e (6, <5i) and t e (t*,oo) such that \u(t)\ = 6 > 6. But \u{t*)\ = 6 so by the 
continuity of u(t), for each 6 e (6,6) there exists t e ( t *, t) such that \u(t)\ = 6. 

It follows from Assumption 2.1 and equation (2.1) that r(6) is a continuous function of 6 and 
hence be Lkr(t)) is a continuous function of 6 as well. Thus we can choose 6 e (6, hi) sufficiently 
small so that A6e ~ Lkr(<5) < 6e~ Lkr ^ for all 6 e (6, 6). 

For such a 6 , noting that u(t) e C k+1 c C 1 for t ^ f, and \u(t)\ > \u(t*)\ = 6, there must exist 
t** e ( t* ,t) such that ^\u(t**)\ > 0 and \u(t**)\ — 6** = sup^^ \u{t)\ e (6,6). This solution 
escapes the ball B( 0, 6**) for the first time at f* > t* > to + kr(6). Let x - u(t**) so |x| = 6** and 

°< iiMn\ 2 = u(n ■ u(n, 

= u(t**) ■ f(f\ u(f*), u(t** - Ti (t*\ u(f*))),u(f * - r,v(r, u(f*)))), 

= x ■ f(f*,x, u(f* - T\(t**, u(f*))),u(f * - T N (f*, u(t**)))). (2.7) 

Since |^(^)| ^ A6e ~ Lkr(lS) < §** e - Lk A5”), follows from Lemma 2.3 part II that f* ^ to+kr(6**). 
Now consider the auxiliary ODE problem (2.4) for I e {1, ..., N} such that either (2.5) or (2.6) 
holds. Let v(9) - u(t** + 9), and noting that [-r/(f**, x), 0] c [-r(h**),0], let 77 ,(0) = u(t** + 9 - 
r i(t** + 9, u(f* + 0))) for 9 e [-r(cT*), 0]. Then 

f* +9- Ti (f* + 9, u(f* + 9)) ^ t 0 + kr(6 **) - 2r(5**) = t 0 + (k- 2 )r(6**), for 9 e \-r(b"), 0], 
It follows that rji e C* _1 ([-r(<T*), 0], B( 0, 6**)) for i = 1,..., N. From (2.7) we deduce that 

x-f(r,x,T 1 m,...,T lN m>o. ( 2 . 8 ) 

Furthermore, we can consider 

77,(60 = U(t 0 + kr(6**) + 9- r,(f 0 + k r (6"*) + 9, U(t 0 + kr(6 **) + 9)), for 9 e [-7-(<f*),0], 
where U is the solution to (1.1) with initial function <t> defined as 


<F(. S ) = 


1 p(s - kr(6**) + f*), s^t 0 + kr(6**) - f*, 
u(s - kr(6 **) + t**), 5 6 [t 0 + kr(6 **) - f*, t 0 ]. 











Thus, (771,..., t]n) g E(k){ 6 **, x, /**). Moreover, 

v(0) = M(f * + 0), 

= /(r + 0, m(T* + 0), «(f + 0 - n(T* + 0, u(t** + 0))),.. 

= /(r +0,v(0),7 7l (0),...,7 7 , v (0)). 


, M(f* + 0 - TA,(r + 0, u(f* + 0)))), 


Also, 

v(-T,(f\ X)) = u(f * - T/(f*, x)) = M(f* - T[(t** , t((f ))) = 77/(0). 
Thus v is a solution to the ODE system (2.4) with / = / and (771,..., 77, v ) 6 £’,/,)(<5 


v(0) = u{t**) - x with |x| = b** so ^7v(.r, 771,..., ?7iv)(0) ■ x - 6* 


,x,f”). But 
which contradicts (2.5). 
Meanwhile, equation (2.8) implies that v(x, 771 ,..., ?7,v)(0) • x > 0 which contradicts (2.6). Thus 
any solution which escapes the ball 5(0, c)”) violates the conditions of the theorem, so u(t) e 
5(0, h**) c 5(0, 6) for all t ^ to. But this is true for all 6 e (b, b) and hence if |</>(.s)[ < A5e~ Lkr(s ' ) 
then \u(Jq, tp)(t)\ ^ b for all t ^ to, establishing Lyapunov stability. Since this holds for all 
A e (0,1) the result follows. □ 

The proof of Theorem 2.5 is complicated by the auxiliary ODE (2.4) being nonautonomous. 
The solution of a nonautonomous ODE escaping the ball 5(0,b) for the first time at t* neither 
implies that j- t \u{t*)\ > 0 nor that there exists f* > f such that \u(f)\ > 5 for all t e (f*, t**). As 
an illustration of this, consider the function y(t) = 6 + t 3 si n 2 (2 tt/ 7) which is easily seen to be 
continuously differentiable and crosses 6 at t = 0 with y'(0) = 0, and for which there does not 
exist any s > 0 such that y(t) > 6 for all t e (0, s). 

Our second main result is to show asymptotic stability of the steady state u - 0 if the auxiliary 
ODE (2.4) satisfies the strict inequality (2.6). We will do this for autonomous DDEs, and for 
simplicity of notation we only present the derivation for problems with one delay term (N = 1). 
The extension to multiple delays is straightforward, and we discuss the extension to periodically 
forced nonautonomous DDEs after Theorem 2.7. Hence we consider autonomous DDEs of the 
form (1.1) with N - 1 for which f(t, u, v) = /(0, u, v) and T\(t, u ) = r(0, u) for all t. In this case 
we may set to = 0 and rewrite ( 1.1 ) as 

u(t) - /( 0, u(t), u{t - t(0, k(1)))), t ^ 0, 
u(t) = (fi(t ) e C, t ^ 0, 


(2.9) 


( 2 . 10 ) 


where C - C([-r(cV), 0], 1R' 7 ). By Assumption 2.1 item 2 the DDE (2.9) has the trivial steady 
state solution u - 0. We write the solution to (2.9) as u(ip)(t) when we want to emphasize 
initial conditions, or just u(t) otherwise. For equation (2.9) the auxiliary ODE introduced in (2.4) 
becomes 

v( 0 ) = /( 0 , v( 0 ), 77 ( 0 )), 0 6 [—r( 0 , x), 0 ], 
v(—r(0, x)) = 77 ( 0 ), 

and since we consider a single delay, there is one such auxiliary ODE associated with (2.9). We 
write the solution of ( 2 . 10 ) as v(x, ? 7 )( 0 ) if we want to emphasize the dependence on x and 77 , 
or just v(0) otherwise. The sets x, t ) defined in (2.3) are no longer dependent on t for an 

autonomous DDE, and so we denote them by E^(6, x) which for (2.9) is defined as follows. 

Definition 2.6. Suppose that Assumption 2.1 is satisfied for (2.9) and k ^ 1. Let 6 e (0, ho] and 
|x[ = 6. Define the set 

( 77 : 77 e C* - 1 ([-t(0, x), 0], 5(0, h)), such that x ■ /(0, x, 77 ( 0 )) ^ 0, j 

E(k)(6, x) = < and for some initial function if e C the solution u(t) of (2.9) satisfies >. (2.11) 
( ?7(0) = u(kr(6 ) + 0 - r(0, u(kr(5) + 0)) for 0 6 [-t(0, x), 0] j 
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To show asymptotic stability of the zero solution to (2.9) it is sufficient to strengthen the 
conditions of Theorem 2.5 by requiring that solutions of the auxiliary ODE problem satisfy the 
strict inequality (2.5) and not the weaker condition (2.6). For the DDE (2.9) the condition (2.5) 
becomes 

—v(x, t/)(0) ■ x < 5, (2.12) 

d 

and we now show asymptotic stability when all solutions of the auxiliary ODE satisfy (2.12). 

Theorem 2.7 (Asymptotic stability). Suppose that Assumption 2.1 is satisfied for (2.9). For 
5 e (0, <5o], x 6 \x\ = 6, define E( k f6,x) as in Definition 2.6. If there exists di e (0, do] such 

that for all 6 e (0, di), and for every x such that \x\ — 6, for all ?/ e E( k )(d, x) the solution v(x, rj)(0) 
of the auxiliary ODE problem (2.10) satisfies (2.12) then the results of Theorem 2.5 hold and 
moreover, the zero solution of (2.9) is asymptotically stable. Furthermore, if\ip(s)\ < 6\e^ Lkr{ ' 6 ' > 
for s 6 [—r(di), 0] then u(t) —> 0 as t —> oo. 

Proof. The only differences between the conditions of Theorem 2.5 and Theorem 2.7 is that 
Theorem 2.5 allows a finite number of delays and nonautonomous / and requires the solution of 
the auxiliary ODE problem (2.4) to satisfy (2.5) or (2.6), while Theorem 2.7 assumes autonomous 
/, one delay, and requires that the strict inequality (2.12) hold. Thus it trivially follows that the 
requirements of Theorem 2.5 are satisfied, and the results of Theorem 2.5 hold. 

Let 5 e (0, di), r = r{§) and |^(i')| ^ 6e~ ,kr for s e [—r, 0]. Then by Theorem 2.5 we have 
\u(t)\ ^ d for all t 0. Consider such a solution. Since \u(t)\ f 6 for all t 0 we have 
lim sup,. <0O \u(t)\ = 6oo with d<x> e [0, d], and it remains only to show that do, = 0. 

Since lim sup, \u{t)\ = doo and {u : \u\ - doo] is compact in there exists t ,■ such that 
lim,--.oo h = oo and lim,- uit,) = Xoo with |x ro | = doo- Assume without loss of generality that 
ti (k + 1 )r for all i. 

Since \u(t)\ f 6 for all t 0, and |</>(f)| 5 ) Se Lkr ^ d for t f 0 it follows from Assumption 2.1 
items 2 and 4 that |Jjw(r)| ^ (To + Ti)d for all t ^ 0. 

Now, consider the sequence of functions vfff) = u(tj + 6) for 9 e [-(k+ l)r, 0]. These functions 
and their derivatives are uniformly bounded with || Vi 11 < 6 and IIJjVj|| ^ (Z^o + L\)d. The set of 
all C 1 functions satisfying these bounds forms a unifromly bounded and equicontinuous closed 
family of functions defined on compact set [-(A: + 1 )r, 0]. By the Arzela-Ascoli theorem the 
sequence of functions v,(0) has a uniformly convergent subsequence. Let {v,| now denote this 
subsequence and let v(6) be the limiting function, which has v(0) = Note that |v(0)| ^ doo 
for all 9 e [-(k + l)r, 0], since the existence of a point with |v(6)| > doo would contradict that 
lim sup,^oo \u(t)\ = doo- 

Let ip*{9) = v(-kr + 9) for 9 e [—r, 0] then we claim that the solution of (2.9) with initial 
function ip* is u*{t ) = v(t-kr) for t e [0, kr]. To see that this is true, let sup, e[0 kri \u*(t)-v(t-kr)\ = 
e ^ 0. Now let uft) solve (2.9) with corresponding initial functions <pj(9) = vf-kr + 9) for 
9 e [-r, 0], so aft) = vft - kr) for t e [0, kr]. For all i sufficiently large we have sup f£ [ 0 kr j \uft) - 
v(t - kr) | = sup fe [ 0i:r ] |vft - kr) - v(f - kr)\ ^ |e by the uniform convergence of the v,- to 
v. But also by the uniform convergence for all i sufficiently large we have \ipf9) - <p*(9)\ - 
\v,(—kr + 9) - v(-kr + 0)| ^ \ j se^ Lkr for all 9 e [-r, 0], and hence by Lemma 2.3 part I we have 
su P/e[(Ur] l«/(0 - M *(0I = sup, e[(Ur] |vft - kr) - u*(t)\ ^ ie. But now 


e - sup | u{t) - v(t - kr)\ < sup |v,-(f - kr) - u*(t)\ + sup [v,(f - kr) - v(t - kr)\ = |e, 

t€\0,kr] te[0,ir] t£{0,kr\ 
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which can only be true if e — 0 so the solution of (2.9) with <p*(0) = v(—kr + 9) for 6 e [—r, 0] is 
indeed u*(t) = v(t - kr) for t e [ 0 , kr\. 

Now let t]( 9) = v{6 - r(0, v(0))) for 9 e [—t(0, Xoo), 0] which implies that r/(9) = u*(kr + 9 - 
T(0,u*(kr + 9))). Moreover |v(0)| ^ 6oo for all 9 e [-(A: + l)r, 0] implies that |//(0)| ^ 6oo for 
9 e [-t( 0, Jfoo), 0] and hence r\ e C* -I ([-r(0, atoo), 0], 5(0, 6oo))- To show that tj 6 £(*)(<5<x>, Xoo) it 
remains only to show that x<x, • /(0, x M , rj( 0)) ^ 0. But if this is false then 

0 > Xoo • /(0, Xoo, 77 ( 0 )) = v(0) • /(0, v(0), 77 ( 0 )) = u(kr) • /(0, u(kr), u\kr - r( 0, u(kr)))) 

= u(kr) ■ u(kr) = \±\u*(kr)\. 

But, \u*(kr)\ = 6oo and j f \u*(kr)\ < 0 implies that there exists s > 0 such that \u*(t)\ > f° r 
t e (kr - e, kr), or equivalently |v(f)[ > for t e (-e, 0). But this contradicts |v(0)| ^ 6^ for all 
9 e [-(k + l)r, 0], so we must have Xoo • /(0, Xoo, J](0)) ^ 0 and 77 6 £(*)( 60c, Xoo). 

Now v(0) = Xoo implies v(0) • x cx , = 6^. But unless = 0 this contradicts that (2. 12) holds for 
all 6 e (0,5i). The result follows. □ 

Notice that Theorem 2.7 not only establishes asymptotic stability of the steady state, but also 
shows that the basin of attraction of the steady state contains the ball 

w : IMI < (2.13) 

We will consider the basin of attraction of the steady state of the model problem ( 1 .2) in Section 6 . 

The extension of Theorem 2.7 to multiple delays is straightforward. The proof given above 
would not be valid for nonautonomous DDEs. However the proof would only fail in one crucial 
step; for a general nonautonomous DDE (1.1), the limiting function v(t) would not in general de¬ 
fine a solution of the DDE. The result is easily extended to periodically nonautonomous DDEs by 
choosing the initial sequence f, to be t, = (k+ 1 )r+iT where T is the period of the nonautonomous 
function /, and if necessary taking a subsequence so that u(t;) converges to x„. 

Our asymptotic stability result and its proof difFers very significantly from other asymptotic 
stability results for RFDEs which are all similar to Theorem 4.2 of Hale and Verduyn Lunel 
[15]. Beyond the technical differences in continuity assumptions, and whether delays are locally 
or globally bounded, there are two fundamental but related differences between our result and 
results such as those in [15]. Firstly, in Theorem 2.7 we establish asymptotic stability, but in The¬ 
orem 4.2 of [15] the stronger property of uniform asymptotic stability is obtained. But secondly, 
auxiliary functions with specific properties are required (in Theorem 4.2 of [15] four auxiliary 
functions, u, v, co and p appear) to obtain the contraction that leads to the uniform asymptotic 
stability. Construction of such functions is difficult even for constant delay DDEs, and a major 
obstacle to the application of these theorems. In contrast, we use a proof by contradiction which 
shows that there does not exist a trajectory which is not asymptotic to the steady state. The con¬ 
tradiction argument establishes asymptotic stability rather than uniform asymptotic stability, but 
does not require any troublesome auxiliary functions, and thus is much easier to apply. In the 
following sections we will use Theorem 2.7 to study the asymptotic stability of the steady state 
of the model state-dependent DDE (1.2). 

We next define the larger sets containing E^)(6, x ) in which we will later show that conditions 
of Theorem 2.7 hold to establish asymptotic stability for the model problem ( 1 .2). By items 4-6 
in Assumption 2.1, if a bound on u(t) is given for t e [-rfhj, ( k - I )r(c))] we can also find bounds 
on up to the k— 1 order derivatives of u(t- r(0, u(t)) for t e [(£— 1 )r(h), kr(6)\. These bounds can 
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be derived from the bounds on /, r and their derivatives. Recalling the definition of E(k)(S, x) in 
Definition 2.6 this leads us to the following definition. 


Definition 2.8. Suppose that Assumption 2.1 is satisfied for (2.9) and k ^ 1. Let h e (0, ho] 
while \x\ = h. Let the functions l)j(d) be Lipschitz continuous in h for j = 0,..., k— 1 and satisfy 


2Dj(<5) > sup 

fE[(A-l)7t<SU-r(<5)] 


u{t - r(f, k(0))|, 


(2.14) 


given that |n(f)l ^ h for all t e [—r(<5), kr(6)\, where u(t) is a solution to (2.9). Define the set 


£(/t)(h, x) 


rt : 77 e PC k -\[-T( 0, x), 0], B(0, h)), x ■ /(0, x, 77 ( 0 )) ^ 0, 
\£m\ < £>j(S) for e e [ r(0, x), 0]J = 0, 1 


(2.15) 


'ik— 1 


where PC k 1 ([—r(0, x), 0], B( 0, £)) denotes the space of C k 2 functions which are piecewise C 

Clearly, E^)(6, x) c £^)(h, x). It is convenient to consider piecewise C k ~ 1 functions in Defini¬ 
tion 2.8 because we will later seek the supremum of an integral over the set S(k)(6, x). Even if all 
the functions in &(k)(6, x) were C k ~ 1 , in general the maximiser could still be piecewise C k 1 . 

In Section 4 we derive bounds Dj{6) for the model problem (1.2), and use these to identify 
parameter regions for which all rj e £^)(h, x) satisfy ( 2 . 1 2 ), and hence the steady state of ( 1 . 2 ) is 
asymptotically stable by Theorem 2.7. For 6 e (0, ho], xeR' f and |x| = h, it is useful to define 


1 

@{5, x) = sup -v(x, t))(0) ■ x, 

i)ES( k) (6,x) d 


T(h) = sup @(6,x), 

\x\=S 


(2.16) 


where v(x, 77 ) is the solution to (2.10). Notice that for h e (0, ho] and |x| = h we have 
sup - v(x, 77 XO) • x < sup - v(x, 77 XO) • x = Q{5, x) < sup Q(6, x)-T (h). 


(2.17) 


rj€E m (S,x) 1 


7lES (k) (S,x) 1 


hi =6 


Thus if T(S) < h for all h e (0,hi) then (2.12) holds for all h 6 (0,hi) and Theorems 2.5 and 2.7 
can be applied. Although T{5) < 6 is a somewhat stronger condition than (2.12) we will find it 
convenient to work with when considering the model problem ( 1 . 2 ). 

The set E ( k)(S,x,t ) given by (2.3) for the DDE (1.1) can be easily generalised to a larger set 
fi(*)(h, x, t), in a similar manner. For t ^ fo + kr(6) we let 


£(/t)(h, x, t) = ■ 


(/,!, ...,i 1n ):t ][ € PC k ~ l ([-r(6), 0], B{ 0, h)), 
x ■ f(t, X, 77! (0),..., ?7;v(0)) > 0, 

\wM0)\ < D, ; (h, r) for 6 e [-r(h), 0], i = 1, •.., N, j = 0,..., k - 1 


where for all solutions u to ( 1 . 1 ) which satisfy |m(s)| ^ h for s e [t - (k + l)r(h), r], 

Dij(6, t ) ^ sup \$ju(s - Tt(s, k(s)))|. 


(2.18) 


(2.19) 


It follows that E(k)(6, x, t ) c £^)(h, x, t), and hence establishing properties on the set &(k)(5, x, t ) 
is sufficient to apply Theorem 2.5. However, we will consider the autonomous model problem 
( 1 . 2 ) in the following sections, and so will not need to consider E(k){6, x, 1) or C ( / f) (h, x, t) further. 
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3. Model Equation Properties 


In the following sections we will apply the Lyapunov-Razumikhin theory of Section 2 to the 
model state-dependent DDE given in (1.2). In this section we consider the properties of the DDE 
(1.2) and its auxiliary ODE (2.10), and will define the sets and functions that we will use to 
apply our results to this model problem. We begin by considering boundedness and, existence 
and uniqueness of solutions of the DDE ( 1 .2) with yU + cr < 0, which generalise the results of 
Mallet-Paret and Nussbaum in [35] for cr < p < 0. 

Lemma 3.1. With c + 0, let p + cr < 0 < a and suppose u 6 C*([0, oo), R) solves (1.2) for t f 0 
with ap( 0) ^ —a then t — a — cu(t) < t for all t > 0. 

Proof. The model DDE (1.2) is invariant under the transformation u i-> —u, c i-> -c, so we 
consider only the case c > 0. Suppose ^j(0) > -a/c and let f > 0 be the first time for which 
u(t*) = -a/c. Then u(t ) > -a/c for t < f implies u(f) ^ 0, but from (1.2) with u(f) = -a/c we 
have ii(f) — (p + cr)u(f) = —~(p + cr) > 0, supplying the required contradiction. If ^(0) = -a/c 
then i/(0) > 0 and the result follows similarly. □ 

We will always consider the DDE (1.2) with a > 0 and p + cr < 0, then Lemma 3.1 assures that 
the deviating argument is always a delay. The lemma also gives the lower bound u(t) > -a/c on 
solutions when c > 0 (or an upper bound on solutions when c < 0). When p < 0 we can bound 
solutions above and below. It is convenient to define 

a acr 

Mq = —, Nq - —, t - a + cN, to - a + cNq. (3.1) 

c cp 

We will use [Mo,lVo] and also [M,N] as bounds on solutions of the single delay DDE (1.2) (in 
contrast to the multiple delay DDE (1.1) for which we used N to denote the number of delays). 

Lemma 3.2. Let p + cr < 0 < a and p < 0 and suppose u e C^fO, oo),R) solves (1.2) for 
t > 0. If cr > 0 let sign(c)M 6 sign(c)[MQ, 0) and sign(c)N > 0, and suppose that sign(c)<p(t) e 
sign(c)[M, N] for all t 6 [-r, 0]. If cr ^ 0 let M = Mq and N = maxjfVo, </?(0)} and suppose 
sign(c)<p(t) ^ sign(c)M{) for all t e [—r, 0]. Then 

sign(c)u(t) e sign(c)(M,N), Vf > 0. (3.2) 

Proof. Again, we consider the c > 0 case, then it is sufficient to show that u(t) > 0 if u(t) - M , 
and ii(t) < 0 if u(t) — N given that u(s) e (M, N) for s 6 (0, t). The case where u(t) — Mq is dealt 
with in the proof of Lemma 3.1, the other cases are straightforward. □ 

Theorem 3.3. Let p + cr < 0 < a. Let the initial history function ip be continuous and for 
p < 0 satisfy the bounds given in Lemma 3.2. For p ^ 0 let sign(c)ip(t) ^ sign(c)Mo for all 
t e (— oo, 0]. Then there exists at least one solution u e C*([0, oo), R) which solves (1.2) for all 
t ^ 0. If p < 0 any solution satisfies the bounds (3.2), while if p ^ 0 any solution satisfies 
sign(c)u(t) ^ sign(c)Mo for all t Q. If ip is locally Lipschitz the solution is unique. 

Proof. Local existence and uniqueness follows directly from the results of Driver [8], and for 
p < 0 global existence and uniqueness follows from the extended existence result of Driver [8] 
using the bounds on the delay and solution given by Lemma 3.1 and 3.2. The only delicate case 
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delay-dependent 


<T 

Figure 2: The analytic stability region £* in the (ji, cr) plane, divided into the delay-independent cone and the 
delay-dependent wedge £ w and cusp £ c . 


is for -cr > /j. > 0 for which (considering the case c > 0) Lemma 3.1 gives only a lower bound, 
u(t ) Mg. But then u(t) ^ nu(t) + ctMq and the Gronwall lemma implies that 

u(t) < (^(0) + - f M 0 = (ip(0) - N 0 V + N 0 . (3.3) 

Since tp(0) ^ Mq > Nq in this case, solutions cannot become unbounded in finite time, and 
global existence again follows. For this case ip(t ) should be defined for all t ^ 0 since with the 
exponentially growing bound (3.3) on u(t) it is possible that t - a- cu(t) —> —oo as t — > +oo. □ 

As already mentioned in the introduction, the constant delay DDE known as Hayes equation, 
which corresponds to (1.2) with c — 0 has been much studied. The {p.,cr ) values for which its 
steady state is asymptotically stable when a > 0 and <: - 0 are well known (see eg. [15]) and 
given in Definition 3.4. 

Definition 3.4 (Stability region 2*). Let a > 0 and c — 0. Let £* be the open set of the (//, cr)- 
parameter space between the curves 


= {(.S', -s) : s e (-oo, 1/a]}, g* = {(ji(s),cr(s)) : s e (0,7r/a)} 

where the functions /u(s) and cr(s ) are given by 

H(s) = scot(as), cr(s) = -■scsc(a.s). (3.4) 

The stability region £* is further divided into three subregions: the cone X.\ = {( 41 , cr) : |cr[ < -//). 
the wedge X w = (2* \ 2a) n {/j. < 0} and the cusp 2 C = 2* n {/j ^ 0], which are shown in Figure 2. 

2* is the parameter region in the ( 41 , cr)-plane for which the zero solution to the DDE (1.2) is 
locally asymptotically stable in both the constant and state-dependent delay cases. The cone 2a 
forms the delay-independent stability region (because this does not change when a is changed) 
while 2 w U 2 C is often referred to as the delay-dependent stability region. For the constant delay 
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case (c = 0) this region is found from the characteristic equation [9]. The results of Gyori 
and Hartung [12] show the state-dependent case (c + 0) of (1.2) has the same (exponentially) 
asymptotic stability region. On the boundary of E* the steady-state is Lyapunov stable for the 
constant delay case, and the stability is delicate in the state-dependent case [40]. 

In this paper we derive new proofs of stability in parts of E* for the state-dependent case 
using Theorem 2.7. The asymptotic stability of the zero solution to (1.2) in all of the delay- 
independent region ((/c, cr) e Ea) will be shown in Theorem 4.1. In Theorem 4.7 we will also 
show asymptotic stability of the steady state of the model problem (1.2) for (/c, cr) in subsets of 
E u UEc by applying Theorem 2.7 with k - 1 to 3. Here we define some notation that will be 
required. Let (//, cr) e E, v U E c , k e Z, k ^ 1, <5o e (0, |a/c|) and 5 e (0, do). It is easy to see that 
Assumption 2.1 is satisfied for (1.2) with Lq = [u|, L\ = |cr|, r max = a and r{8) — a + |c|d. Thus 
for the model problem (1.2) the sets E^( 6 , x) from Definition 2.6 are given by 

( 77 : rj e C k ~ ] ([-a - cx, 0], [-d, 5]), fj.x 2 + crxr](0) ^ 0, j 

E(k)(d, x) — < and for some initial function if e C equation (2.9) has solution >. (3.5) 

( 77 ( 0 ) = u(kr(jS) + 6 - a - cu(kr(5) + 9 )) for 9 e [-a - cx, 0] J 


To apply the stability theorems in the next section we will derive bounds Dj(6) for j = 0, 1 

and 6 6 (0, <5o] as in Definition 2.8. Once these bounds are determined, the sets S^)(S,x) from 
Definition 2.8 are given by 


S(k)(8, x) = 


77 : 77 6 PC k '([-a - cx, 0 ], [- 5 , < 5 ]), fix 2 + crxr]( 0 ) ^ 0 , 1 

|^j77(0)| ^ £>j{6) for 9 e [-a - cx,0], j = 0,..., k - 1 J ' 


Let r+ — a + c 6 then the auxiliary ODE problem (2.10) becomes 


v{9) = iuv(9) + cnj(9), 9 e [-/+, 0], 

v(-r+) = 77 ( 0 ). 


(3.6) 


(3-7) 


Integrating (3.7) yields. 


r° 

v(0) = 77(0)e / ' r+ + cr I e^ e r 1 {9)d9. 


(3.8) 


Since the DDE (1.2) is scalar the set of x such that [x| = 6 consists of just two points x - 6 and 
x — - 6 . Suppose first that x — 6 , then (3.6) implies that 77 ( 0 ) e [-<5, -dfj /<r]. 


Definition 3.5. Let a > 0, c^0, cr^/r and cr < - 77 . For any 6 6 (0, |a/c[) and u e [-<5, -6/u/cr], 
define r+ = a + c6 and rj(k)(9 ) for 9 e [-r+, 0] by 


9 (k) ( 8 ) = jnf 77 ( 0 ). (3.9) 

(W) 

!?( 0 )=« 

We also define the function I(u, 6 , c, k) to be 

r° 

I(u, 6 , c, k ) = U(T+ + cr e^ e r 1 (k) {9)d9. (3.10) 

J-r, 

The function 77 ^) given by (3.9) is the most negative one in d) satisfying 77 ( 0 ) = u, and so 
since cr < 0, this function maximizes v(0) for fixed 77 ( 0 ) by maximising the second term in (3.8). 
This is the reason for considering 77 e PC k ~'([-u - cx, 0], [-5, c>]) in the definition of £>( k ){5, x). 
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We really want to maximise v(0) for ?/ e E^fd, x), where the smaller set E^(6, x ) is defined in 
(2.11). Even though all the functions 77 6 E ( k)(6,x) satisfy 77 e C k ~ l ([-a - cx, 0], [-5,5]), the 
maximiser will in general only be piecewise C* . With Definition 3.5 we can derive bounds on 
the solution v(0) of the auxiliary ODE (3.7) for all 77 e £>(k){6, x) in both cases where x — ±5. 

Lemma 3.6. Let a > 0, c 4- 0, <r F jj and cr < —jj. Let 6 e (0, |a/c|). The solution of the auxiliary 
ODE system (3.7) satisfies 

v(0) < sup J(w, 5, c, k), V 77 e £(*)(5,5), v(0) > - sup I(u,6,-c,k), V 77 e £^)(5,-5). 

Proof. First consider x — 5. The function I(u,6,c,k ) comes from (3.8) and depends on c and 
5 through r+. Since, as noted above, the choice of 77 ^) maximizes (3.10) for fixed u, the first 
inequality in the statement of the lemma follows. 

Next consider x — -5, then (3.6) implies that //CO) e [5/t/cr, 5]. This time we should consider 
the most positive function in & (^(5, -5) satisfying rj(0) - it e [6/i/cr, 5], to obtain a lower bound 
on v(0) for all 77 6 (5, -5). However, the model DDE (1.2) is invariant under the transforma¬ 
tion (m, c) h-> (-«, -c), so this function is and the second inequality follows. □ 

Notice from (3.10) that the functions I(u,6,c,k) and I(u,5,-c,k ) only differ in their inte¬ 
gration limits with I(u,6,c,k) integrating 77 ^) over the interval [-a - c5,0] and I(u,6,-c,k) 
integrating over [-a + c6, 0]. The integration over the larger of these intervals will be important 
in the following sections and so it is convenient to define 

P(6,c,k)~ sup I(u,6,\c\,k). (3.11) 

ue[-6,-6ji/<r] 

Comparing the cases when x = 5 and -5 has to be done separately for each value of k. and we 
can also explicitly define the functions rj^ for each k. This is handled in the following section 
where we show that P(5, c,k) < 5 implies 'F(d) < 5, and apply Theorem 2.7 to obtain asymptotic 
stability for {(//, cr) : P(5, c, k) < 5}. 

Barnea [1] applied Lyapunov-Razumikhin techniques to the c - 0 constant delay case of the 
model DDE (1.2). His results do not apply to state-dependent case, as they were based on a result 
for autonomous RFDEs which assumed F was Lipschitz, and he did not define an auxiliary ODE, 
nor sets similar to E^)(6,x) or <3,^(5, x). However, he did define functions i],^ for the constant 
delay case by considering the most negative bounded functions with k — 1 bounded derivatives 
as the function segments in the RFDE. In the limit as c —> 0 our ?/, 7 0 functions reduce to those 
found by Barnea for the constant delay case. Because of the linearity of ( 1 .2) with c = 0, Barnea 
did not have to consider the upper and lower bounds separately as we did in Lemma 3.6, but 
did define a function which is equivalent to P(5,0, k) in (3.11). Our asymptotic stability results 
for the state-dependent model DDE (1.2) constitute a significant generalisation of the Lyapunov 
stability results of Barnea [1] for the constant delay case, and moreover in Section 5 we will 
correct an error of Barnea for the k — 2 constant delay case. 

4. Asymptotic stability for ii(t) = fiu(t) + cru{t - a - cu(t)) using E>^)(6,x) 

In this section we consider the model state-dependent DDE (1.2) and use Theorem 2.7 to show 
that the steady state is asymptotically stable in various parameter sets. In Theorem 4.1 we use 
the set £(q(5, x) to show that the steady state is asymptotically stable whenever the parameters 
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values (fi, cr) are in the cone X,\. In the rest of the section we consider parameters in the wedge 
and the cusp (2 W U2 C ), and use Theorem 2.7 with k — 1, 2 and 3 to show the steady state is 
asymptotically stable for (ju,cr) e \P(\, 0. k) < 1), where P( 6 ,c,k ) is defined in (3.11). The sets 
{P(l, 0, k) < 1} are nested in 2 W U 2 C , becoming larger with k. We also find lower bounds on the 
basin of attraction of the steady state. For the constant delay case (c = 0) the parameter regions 
found in U 2 e are independent of the choice of 6 in x). For the state-dependent case, 

these regions change with c and 6 (see Figure 4) and converge to the region for the constant delay 
case as 5 —> 0. 

We begin by showing asymptotic stability in the cone 2 a- The following result could also be 
shown by adapting a stability result for time-dependent delays, such as that of Yorke [43]. Recall 
that Mo is defined by (3.1). 

Theorem 4.1 (Asymptotic stability for (1.2) in Ea)- Let a > 0, c 4 Qand\cr\ < —ji so (ji,cr) e Ea- 

If\ip(t)\ ^ \Mo\for t e [—a — c|Mo|, 0] then the solution u(t ) to (1.2) satisfies u(t) —> 0 as t —> oo. 

Proof. With |x| = 6 and ft < —|cr| it is impossible to satisfy fix 2 + crxrjiQ) ^ 0 with \rj(0)\ ^ 6 
and so £(i)(d, x) is empty and asymptotic stability of the steady state follows from Theorem 2.7. 
This holds for all 6 e (0, |Mq|] and it follows directly from Theorem 2.7 that u(t) —> 0 as t —> oo 
provided |^(f)[ < \Mf\e^ Lr< ^ M °^ for t 6 [-a - c|Mo|, 0]. However, the exponential correction term 
e -Lu\M a \) comes f rom using Lemma 2.3 in the proof of Theorem 2.7 to ensure that \u(fi)\ < \Mq\ 
for t e [0, r(|Mo|)]. But Lemma 3.2 already ensures that |n(t)| < \Mq\ for all r > 0 if |^(f)l ^ |Mo| 
for the model DDE (1.2); the result follows. □ 

We already noted in Section 3 that Assumption 2.1 is satisfied for (1.2), and derived an ex¬ 
pression for£(i)(d, x) and indicated which are the most relevant functions in these sets. For k - 1 
we do not need any bounds Dj( 6 ) and the members of the set Bafd, x) need not be continuous. 
Then the function //, |, from (3.9) is given by. 



(4.1) 


For k - 2 we need to find D\( 6 ) such that 1^77(601 ^ Di(£) for all rj e E ( 2 ){ 6 ,x), where 
E( 2 )(S, x) is defined by (3.5). For;; e E^ifi, x) we have 


tj( 0 ) = u( 2 r + 8 - a - cu( 2 r + 8 )) for 8 6 [—a - cx, 0] £ [-r, 0], 


where \u(t)\ ^ 6 for t e [—r, 2r] and solves (1.2) for r ^ 0. We easily derive that |n(r)| ^ 
\fiu(t)\ + | cru(t - a - cu(t))\ ^ Qfi\ + |cr|)<5 for t 6 [0, 2 r]. Then 

t]'( 8 ) = ^gii( 2 r + 8 - a - cu( 2 r + 8 )) — (1 - cii( 2 r + 8 ))ii( 2 r + 8 - a - cu( 2 r + 8 )). 

Hence \r/'( 8 )\ ^ D\5 for 8 e [—r, 0] where 


D\ = (|/i| + [cr|)(l + (\fi\ + |cr[)|c|5). 


(4.2) 


Thus we can choose 'D\ (d) = D \6 (note that D\ also depends on 6 ) to define )(d, x) and we 
obtain that £’( 2 )(d, x) c &( 2 )(< 5 , x). The function is given by (3.9), as 




8 e [-a - c 6 , when < r + - a + c 6 . 
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(4.3) 


where D\ is defined by (4.2). 

For ij e £( 3 )(d, x) c E( 2 j (d, x), the same bound on the first derivative of // applies, and we also 
bound the second derivative as follows. We have 


T](0) = u(3r + 6 - a - cu(3r + 9)) for 6 6 [-a - cx, 0] c [-r, 0], 


where \u(t)\ ^ 6 for 1 e [—r, 3r] and solves (1.2) for r ^ 0. As above we have that \u(t)\ ^ 
(|yu| + |cr|)(5 for t e [0,3 r\. Now noting that t - a - cu(t) e [0,2r] for t e [r, 3r] we have 

|ii(0l = | fju(t) + cru(t -a - cu(t)( 1 - cu{t))\ < (|/i| + |cr[) 2 (l + \o~c\S) 6 , 

for t 6 [r, 3r], Then, since 3r + 0 - a - cu{3r + 9) e [r, 2/-] for 0 6 [—r, 0] it follows that 

\jj"(9)\ = |^«(3r + 9 - a - cu(3r + 0))| 

= |(1 — cu(3r + 9)) 2 ii(3r + 9 - a - cu(3r + 9)) - cu(3r + 9))u(3r + 9 - a - cu(3r + <?))| 

< ( D\ + (M + |cr|) 3 |c|d)(l + \crc\ 6)6 = D 2 5. 

Hence for all r\ e E^( 6 , x) we have \q"{9)\ ^ r D 2 (S) — D 2 6 where 

D 2 = {D\ + (M + |cr[) 3 [c|<5)(l + |crc|d), (4.4) 

and lim^o D 2 = (lim tf _ > o D \) 2 = (|ju| + |cr|) 2 . Taking D\ and D 2 to satisfy (4.2) and (4.4) ensures 
that 6 ( 3 )(5, x) c £( 3 )(d, x). Then the rj^) function from (3.9) for k = 3 can be defined by 


*7(3)(0) - *7(3) (0 + #shift), 0 e [-r+, 0] 

where 

/2(fi + 5)\* 

\ « 2 h j ’ 

- C <50? 

“ + 5 + 257 

’ 


'7(3) (6) - 


- 5 , 

-5 + 


\D 2 9 2 


2 ' 


9 ^ 0 , 

0 e (0, j£). 


^shift — 


(4.5) 


(5£>t 


Ci 6 [-h, —<5 + 557 ], 


« > -(5 + 


<5Z>f 

2ZV 


(4.6) 


Here (9 s hift is a convenient device which allows us to define ij( 2 ,(0) for all values of u by the single 
function 770 ) ( 0 ) with the shift used to obtain the correct value of u. 

The ijiti functions define I(u, 6 , c, k ) via equation (3.10) and P( 6 , c,k ) through equation (3.11). 
For k = 1, using (4.1) we easily evaluate 


, , (--de^ + d-i \-en, F^0, 

P{ 6 , c, 1) = l(-/i 6 /cr, 6 , |c|, 1) = 1 cr n 

[-dcrr, fx — 0. 

For k = 2, from (3.10) and (4.3), if ^ r+ then 


I{u, 6 , c, 2 ) = ueP r+ + cr 


(u + Di69)d9 


r° 

I e-»% 

J-r+ 

juje^ + Z(eP r + - 1 )] + ZDi6[±(eP r + - 1 ) - r + e>"'], /j + 0 , 

F - 0, 
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I /s ~ aD\r\ r- 

I W + 0T+M- 2 ^^’ 


(4.7) 


(4.8) 


(4.9) 










Figure 3: Sample t]^)(^) functions for a = 1, c = 0 and (5=1. The T](k){0) functions are defined in (4.3) and (4.5). 


while if < r+ then we have to split the integral into two parts and 



e~^(-5)d6 + cr f e~ M \u + D^de 

I 5+u 
^ D\ 6 


m) + ^ 6 \r{e ^ D ' 6 


u - cr6r+ + + m) , 


H + 0 , 

P = o. 


(4.10) 


(4.11) 


To determine P(6, c, 2) we perform the integration in I(u, 6, |c|, 2) and find u to maximise this 
function. If |±| ^ r then u e [(rD\ - 1 )6,-6fi/cr\. This is only possible in the region where 
rD\ - 1 ^ —fi/cr. From (4.9) we have 


I(u, 6, [c|, 2) = Ii(u, <5) := 


_ (u[e» r + f(e^ r - 1)] + - 1) - re^], /i + 0, 


u + crrit — 


crDii 2 


6 , 


l-i - 0 . 


(4.12) 


If < r then u has another upper bound u < (rD i - l)d so u e [-d, min{(rZ>i - l)d, -^d}]. 
Since the integration is broken down into two parts in this case we label the expression we derive 
as I 2 , and from (4. 1 1) we have 


I(u, 6, |c|, 2) = I 2 (u, 6) 


{■’(«"- 5 ) + 

\u-cr6r + jj^(6 + u) 2 , 


l)-^ r ], /j4 0 , 

M = 0. 


(4-13) 


The main differences between the expressions for I(u, 6, c, 2) and I ( u , 6 , |c|, 2) are that the former 
involve r+ = a + cd, and the latter use r — a + c|c) as well as being subject to different restrictions 
on the values of u for which they apply. In (4.9), (4.1 1),(4.12) and (4.13) the // - 0 expressions 
equal the /j —> 0 limit of the // 4 0 expressions. Results for /j = 0 thus follow from those for 
j-i 4 0, and so we do not treat these cases separately below. 


Theorem 4.2. Let a > 0, c 4 0, cr ^ and cr < -/a. Let 6 6 (0, \a/c\). IfP(6, c, 2) < 6 then 


if rD\ - 1 < -p/cr, 
I 2 (-6p/cr,6), if rDi - 1 > -///cr, 
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P{6, c, 2) = 


(4-14) 












where I\ is defined by ( 4.1 2) and I 2 is defined by (4. 1 3). 
Proof. See Appendix A. 


□ 


We will not state an explicit expression for P(d, c, 3). When needed, this can determined by 
evaluating (3.11) numerically for k — 3. 

We now prove four lemmas which will be needed for the proof of Theorem 4.7 where we show 
asymptotic stability in the set {(//, cr) : P( 1,0, k) < 1}. 

Lemma 4.3. Let a > 0, c ^ 0, cr fi p and cr < —p. Let 5 e (0, |a/c|) and it e [—<5, -dp/cr] be 
fixed. Then I(it, 6, c, 1) decreases with decreasing r+. 

Proof. For r+ > 0, 

-^—{ue^ r * + cr J' e~ f,e r] ( i)(0)d0j = e Mr *[pu + crrj (l) (-r + )\ = e^ r+ [pu - cr6\ > 0 , 

since cr < -p and fi e [—b, -dp/cr\. □ 

Lemma 4.4. For k = 2 or 3, let a > 0, c + 0, cr ^ p and cr < -p. Let 5 e (0, \a/c\) and 
u e [—<5, —pd/a] be fixed. Let I(r+) be the expression for I (fi, 6, c, k ) as a function of only r+; 

r° 

I(r+) = ue^ + cr I e~^ e r] (k) (6)d6, ^J(r+) = eP r *\pu + arp k) (-r + )\. (4.15) 

J-r+ 

(A) Ifp fi 0, then £'-I(r+ ) > 0. 

(B) Ifp > 0 and rj^ k )(—r + ) fi 0, then -J^-I(r + ) > 0. 

(C) If -£-I(r+) fi 0, then p > 0, ppf-rf) > 0, and I(u, 6, c, k ) < 6. 

Proof. Parts (A) and (B) are easy to show. Let -£-I(r+) fi 0. From the first two cases, this is 
only possible if p > 0 and tjp)(-r + ) > 0 . 

Consider & = 2 first. Since Tj( 2 )(-r+) + -5 we are in the case jjfi > r+. Thus I(u,d,c, 2) is 
given by (4.9). Since r/( 2 )(-r + ) — it - D\dr + > 0 implies — £ D\dr+ < jh, we deduce 

I(u, 5, c, 2 ) < u\<f r + + -(<fi r + - 1)1 + d^(e^ - 1 ) + -uef"*, 

1 p J p~ p 

= -fi— + (1 + 2 —)ue lir * + b^-^-(e / "' + - 1 ), 

p p p- 

fi -fi— 1 — ue IJr+ + d — ir-(ef r * - 1 ), since — < -1 and fi ^ rjn\(—r + ) > 0 , 

P P P- P 

= -(e^ - \)(u + ^-d) <0 <d. 

p » 

For k - 3, from D^/D 2 ^ 1 it follows that 

dD 7 d d 

hO)(Di/D 2 ) - -d + -jjj- fi -d + - - -- < 0. 

Since ?/( 3 ) is an increasing function and we require t](y)(-r+) > 0, then —r+ + 0 s hif t > § 7 . Thus 
9 + 0 s hift > D 1 ID 2 for all 9 e [-r+, 0]. By the definition of 77 ( 3 ), in this case tjq)(9) - tj( 2 ){ 9 ) for 
9 e [-r+, 0] and J(fi, d, c, 3) = J(fi, d, c, 2) < d. □ 
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Lemma 4.5. For k — 1, 2 or 3, let a > 0, c + 0, cr ^ /r ancl cr < -ji. Let 6 e (0, \a/c\). Then 


fi,cr) : P{5,c,k)~ sup I(u,6,\c\,k) < 6\Q l{p.,cr) : sup I{u, 5, -|c|, k) < 6 

ue[-d,-5/j./o-] ) V ue[—6,—6fi/cr] 

Proof. Recall from Section 3 that changing the sign of c in I(u, 6, c, k) only changes the value of 
r+ — a + c6. For k = 1 the result follows from Lemma 4.3. 

For k — 2 or 3, let e {P(6,c,k ) < 5} and u e [—6,—6fi/cr]. Recall that cr < 0, while 
r](k)(Q) is a nondecreasing function in 9. There are two cases to consider: 

(i) If juu + - |c|5)) < 0 then by Lemma 4.4(C), J(w, 6 , -|c|, k ) < 6 . 

(ii) If fiii + o-T](k)(—(a + |c|5)) ^ fiit + crt]^)(.-(a - |c|5)) > 0, then fiu + <Ttj{j.)(-T)) > 0 for all 
r € (a — \c\6,a + |c|<5). By equation (4.15), the expression for I(r + ) is increasing over this 
interval and thus, I(u, 6 , —|c|, k) ^ I(u, 5, |c|, k) ^ P(5, c, k) < 6 . 

Thus I(fi, 6, —|c|, k) < 6 and the result follows. □ 

Lemma 4.6. For k - 2 or 3, let a > Q, c ^ Q, cr ^ ji and cr < -ji. IfO < 5„ ^ < |a/c| then 

{(jj.,ct) : P(6„,c,k) < 6„] c {(jd,cr) : P(6„c,k) < 6*}. 


Proof. Increasing <5 increases r — a + \c \6 which is the only source of nonlinearity in 6 in the first 
expression (4.7) for P( 6 , c, 1). Thus for jj 4- 0 


d (P( 6 ,c, 1) 
551 6 


)= + f > °- < 416 » 


Positivity also follows trivially from (4.7) when /i = 0. The result follows for k = 1. 

For k = 2 or 3, consider I(s 6 , 6 , |c|, k )/6 and note that r, D\, ..., Dk- 1 are the only terms in the 
expression that depend on 6 , and that increasing 6 increases r, D\ and Dt- Thus 


dS\ 


dr l 


k-l 


7=1 


d 

(I(sd, 6 , |c[, k)\ 

90, 

dDj 

l ^ J 

1 d 6 


We focus on the first term on the left-hand side, since all the remaining terms are positive. From 
(3.10) we can write 

_5_ | JU5,5,|f|,<c) j = + = c jir^ s + z^ (jfc) (_( a + | c |5)]. 

Let r* - a + |c|5„, r** - a + \c\5„ and (yu, cr) e { P( 6 „,c,k) < 5,,}- Let s e [-1 ,-fi/cr] and 
use the notation 9) to denote the function t]^) as a function of both 9 and 6 . Note that 
t](k)( 6 , -(a + |c|5))/5 is always decreasing with 6 . Consider the following cases: 

(i) If /us + o-tj(k)( 6 », -r*)/d* ^ 0 then by Lemma 4.4(C), I(s 6 », 6 *, |c|, k) < 6 *. 

(ii) If /js + o~r](k)( 6 t,, —r**)/d** > fis + <TT]^ k )(d t ,-r*)/ 6 , > 0 then ^ 0 for 6 e 

[5„,5„]. Thus, §- s ( I(s 6 f c ’ k) ) > 0 ford 6 [5.,5„] and, 

I(sd,,d„,c,k ) < I(sd„,d„,c,k) P(d M ,c,k ) < ^ 
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-4 0^2 

(a) |P(1,0,1)< 1} 



-4 0 M 2 

(b) {P(1,0,2) < 1| 


Figure 4: For (a) k = 1 and (b) k = 2, the set (P(1,0, k) < 1| is shaded green and the boundary of |P(l/2,1 ,k) < 1/21 
is drawn in red. If (//, it) € {Pi 1.0. k) < 1| then the zero solution to (1.2) is asymptotically stable by Theorem 4.7. As 
5 —> 0 the proof of Theorem 4.7 shows that {P(<5, c, k) < ci} converges to { P( 1.0, kj < 1|. 


Cases (i) and (ii) both yield I(s8*, 8,, |c|, k) < 8,. Since this holds for all s e [-1 ,-ju/cr], 
P(b*, |c|, k) < 8* follows. □ 

With these lemmas we can prove our main result. 

Theorem 4.7 (Asymptotic stability for (1.2) using S(kf 8 , x)). For k — 1, 2 or 3, let a > 0, c A 0, 
and (/u,o~) 6 {P(l, 0, k) < 1} where P{8, c, k ) is defined by (3.11). Then (//, cr) e {P(di, c, k) < hi) 
for some 8 i e (0, |a/c|). Furthermore, for 8 6 (0, hi] let 82 = and \ip(t)\ < 82 for 

all t e [—a — [c|d, 0], then the solution to (l .2) satisfies \u(t)\ ^ 8 for all t ^ 0 and lim^c*, u{t) = 0. 

Proof. For this proof define 

7= (J {P(b, c, k) < d}. 

5e(0,|a/c|) 

First we show that J = {P( 1,0,1) < 1). When c - 0 it is seen that I(s8,8, 0, k)/8 is independent 
of 8 for k - 1, 2 or 3. From this it follows that P(8, 0, k)/8 = P( 1,0, k). Moreover, for all c, when 
8 —> 0 then r —» a, and I(s8,8, \c\,k)/8 —» I{s8,8, 0, k)/8, since c only appears multiplied by 8 
in these expressions. Thus P(8,c,k)/8 —» P(1,0, A:) as 8 0. Because of this and Lemma 4.6, 

J = {P(l, 0, k) < 1). 

Let (ju, cr) e {P( 1,0, k) < 1). The existence of 8\ such that (ju, cr) e {P(8i,c,k) < hi) follows 
from the above discussion. It also follows that (/j, it) e [P(6, c, k) < <5) for all 8 e (0, d) ]. 

Let 8 e (0, hi]. Consider the auxiliary ODE (3.7). For all rj e 8) it follows from 

Lemmas 3.6 and 4.5 that v(0) ^ sup ae ^_ s I(u, 8, e, k) f P(8,c,k) < 8. Similarly, for all 
r\ e S(k)(8, -8) we obtain v(0) ^ - sup ae[ _ a _^ /o .] I(u,8, -c, k) > - P(8,c,k ) > -8. Thus (2.12) 
holds for all rj 6 S(k)(d, x) or any \x\ = 8. This is true for all 8 e (0, hi]. Since E(kf8, x) c £^(5, x), 
applying Theorem 2.7 completes the proof. □ 

For given (//, cr) the condition P(l, 0, k ) < 1 ensures (2.12) is satisfied for all r\ 6 F>(kf8, x) and 
hence Theorem 4.7 establishes asymptotic stability for (ju, cr) in the part of X u . U X r for which 
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P( 1,0, k) < 1. These sets are shown in Figure 4 for k = 1 and k - 2. The stability region 
{P(1,0,1) < 1} shown in Figure 4(a) for the DDE (1.2) comprises a relatively small part of 
Zi, U 2 C , because it is derived by requiring that (2.12) holds for all rj 6 <5, ],( 6 , x). But £(i)(5, x) is 
a very large set, with the main restrictions on rj being that it is merely piecewise continuous with 

INI ^ S. 

We obtain a larger stability region by increasing k. This is seen in Figure 4(b) where 
P(1,0,2) < 1 ensures that (2.12) is satisfied for all // e )( 6 ,x) results in a significantly 
larger stability region than seen in Figure 4(a). Since £ C 2 )(d, x) c £ ( \ ,(<5. a), with all functions 
rj e &( 2 )( 6 ,x) satisfying the derivative bound \jj'( 6 )\ ^ D\5, the set E>( 2 ){ 6 , x) is smaller than 
£(i)(d, x) and it is possible to satisfy (2.12) over a larger region of (//, cr) parameter space. We 
will compare the sizes of the stability regions [P(l, 0, k) < 1) for different k in Section 5. 

The boundary of [P(l/2, 1, k) < 1 /2) is also shown in Figure 4 for k — 1 and k — 2. As 6 —> 0 
the sets {P( 6 ,c,k) < d) converge to the set {P(1,0, k) < 1), and for (j-i,cr ) e {P(1,0, k) < 1} 
the inequality P( 6 , c, k ) < 5 can be used to determine the largest 5\ and hence the largest cL for 
which Theorem 4.7 applies. This determines a ball which is contained in the basin of attraction 
of the steady state, and in Section 6 we consider how the size of this lower bound on the basin of 
attraction varies with k. 

5. Comparison of the stability regions 

In this section we compare the sets in which we can establish asymptotic stability of the steady 
state of the state-dependent DDE (1.2) using Lyapunov-Razumikhin techniques. In Sections 4 
we showed asymptotic stability for (//, cr ) 6 2a, and for (//, cr) in the parts of the cusp 2 C and 
wedge 2 W for which P( 1,0, k) < 1 for k - 1,2,3. Measurements of these sets and the exact 
stability region 2* are presented in Tables 1-3, and they are illustrated in Figure 5. 

To compute these stability regions, from Theorems 4.7 we need to compute P(1,0, k) in the 
limiting case c = 0, 6 = 1. This was done in MATLAB [36]. For k - I and 2 we have exact 
expressions for P(1,0, k) given by (4.7) and (4.14). Noting that cr < 0 in 2 W U 2 C , from (4.7) we 
find that (ji, cr) satisfies P( 1,0,1) < 1 when 

- eT) - <r - vuT > 0. (5.1) 

The boundary of{P(l,0,l)< 1} is defined by equality in (5 .1 ). 

For k — 3 the value of P(l, 0, 3) was calculated by maximizing the function I(u, 1.0, k) over 
u 6 [-1, -fj/cr] using the MATLAB fminbnd function. 

The boundary of [P(l, 0, k) < 1} is then found by fixing one of /j or cr and using the fzero 
function to find the value of the other one which solves P( 1,0, k) - 1 = 0 (except in the case k - 1 
where for given /i, applying the quadratic formula to (5.1) determines cr). The largest value of jj 
for each region (shown in Table 3) is then found by regarding the ji that solves P(l, 0, k) = 1 as 
a function of cr and using fminsearch to find the cr that maximises //. The boundary of the full 
stability domain 2*, found by linearization, is given by Definition 3.4. 

Since by Theorem 4.7, at least for k 3, the Lyapunov-Razumikhin stability regions are given 
by P(l, 0, k) < 1 irrespective of the value of c, we obtain the same regions in the constant c - 0 
and variable c £ 0 delay cases. This is consistent with the linearization theory of Gyori and 
Hartung [12] who showed that 2* is the exponential stability region for both c — 0 and c 4- 0. 

When ji = 0 the DDE (1.2) becomes 


ii(t) — cruit —a — cu(t)) 
23 


(5-2) 



(a) I A U{P(1,0,1) < 1} 


(b) £ A U{P(1,0,2)< 1} 


(c) S A U{P(1,0,3)< 1} 


Figure 5: Stability regions found using Theorems 4.1 and 4.7 shaded in green. Parameter pairs (jJ,cr) e which are 
contained in the wedge or cusp £ c but not in {P(l, 0, k)} are shaded dark or light grey respectively. 


Region 

cr at ^ = —5 

cr dX n = —2 

cr at /r = 0 

[P(1,0,1)< 1} 

-5.0673855 

-2.5578041 

-1 

[P(l,0,2)< 1} 

-5.0676090 

-2.5875409 

-1.5 = -3/2 

[P(l,0,3)< 1) 

-5.0678325 

-2.6127906 

-1.5416667 = -37/24 

2* 

-5.6605586 

-3.0396051 

-1.5707963 = -;r/2 


Table 1: Boundaries of the stability regions: Values of cr for fixed /j. with a = 1. 


and the intervals of cr values in the stability regions, (shown in the last column of Table 1) can 
be found exactly. From (4.7) we have P(l, 0,1) = -era when /j = 0 which implies cr e (— 1/a, 0) 
for P( 1,0,1) < 1. Similarly, when cr < -1/a from (4.14) we have P( 1,0,2) = -era - 1/2 and 
hence cr = -3/(2 a) on the boundary of [P(l, 0,2) < 1). Magpantay [30] shows that for/r = 0 we 
require cr e (-37/(24a), 0) for P(1,0,3) < 1. For the constant delay case of (5.2), with c = 0, 
Barnea [1] showed Lyapunov stability for cr e [-3/(2a), 0], by applying Lyapunov-Razumikhin 
techniques with k = 2. The stability bound cr ^ -37/(24a) seems not to have been derived 
before for the constant delay case of (5.2), but is well-known for Wright’s equation [42] which is 
a nonlinear constant delay DDE whose linear part corresponds to (5.2) with c = 0 (see [28, 42]). 

When /r = 0, the boundary of \P( \, 0, k) < 1} seems to converge rapidly to -n/(2a), the 
boundary of £*, as k — > oo, suggesting that the full stability interval can be recovered. Indeed, 
for constant delay with \i — c - 0 Krisztin [26] showed Lyapunov stability for cr e (—zr/(2a), 0] 
by considering k —> oo. 

Barnea [1] and Myshkis [37] also applied Razumikhin techniques to establish Lyapunov sta¬ 
bility for (1.2) in the case of constant delay (c = 0) with /.i + 0. The regions in which they 
claim stability are shown in Figure 6. The region found by Myshkis [37] has cr ^ -1/a and 
/j ^ cr(acr + 1 )/(acr - 1), which for a = 1 always has /j 5 / 3 - 2 V2 and is contained in 
[P(1,0,1)< 1). 

Barnea [1] claimed that the Lyapunov stability region of (1.2) with c = 0 contains the region 
X 2 = ](>u, cr) : 0 < s* ^ a, P < 1} where 




cr 


P = 


cr(ju + cr) 

P 2 

24 



cr 


s 


yU + cr 


























Region 

fi at cr = — 5 

/r at cr = -2 

/j at cr = — 1 

{R(1,0,1)< 1} 

-4.9286634 

-1.16401463 

0 

{£(1,0,2) < 1} 

-4.9283941 

-1.0040856 

0.38774807 

{£(1,0,3) < 1} 

-4.9281247 

-0.93607885 

0.39925645 

I* 

-4.2734224 

-0.63804505 

1 


Table 2: Boundaries of the stability regions: Values of ji for fixed cr with a = 1. 


Region 

Supremum of y 1 

Corresponding value of cr 

{£(1,0,1) < 1} 

0.18822641 = ln((l + V2)/2)) 

-0.45439453 

{£(1,0,2) < 1} 

0.45697166 

-0.73935547 

{£(1,0,3) < 1} 

0.45700462 

-0.74059482 

2* 

1 

-1 


Table 3: The values of /j and cr at the rightmost boundary point of each stability region with a = 1. 


We show this region in Figure 6(a), but Barnea did not actually graph AN or give its derivation 
in [1]. He noted that setting P — 1 and letting /r —> 0 yields that the point cr — -3/2a is a 
boundary of X 2 on the cr-axis. We observe that setting s* — 0 and /r —> 0 yields cr = -1 /a as 
the other boundary of X 2 on the cr-axis. Thus the region AN does not include the whole interval 
cr e (-3/2a, 0] on the cr-axis which Barnea had proven to be Lyapunov stable in the /j = 0 
case in the same paper [1], Barnea’s stability region AN is hence incomplete. Although the ija) 
function used by Barnea to show Lyapunov stability corresponds to (4.3) with c = 0, it appears 
that Barnea performed his integration assuming that ^ r + in all cases. The case when 
> r + occurs in the n — c — 0 case (as well as the general case c + 0, fi + 0 considered in 
(4. 1 2),(4. 1 3)). Omitting this case results in the incorrect stability region AN- The correct region 
is (ju, cr) e {£(1, 0,2) < 1) as illustrated in Figure 5(b). Moreover within this region we show the 
stronger property of asymptotic stability for both the constant delay (c = 0) and state-dependent 
delay (c + 0) cases. 

Tables 1 and 2 also show that for // < 0 we can show asymptotic stability in a larger part 
of £* by increasing k. However when n 0 the improvement in going from k - 1 to 2 to 3 
is very marginal and we can only show asymptotic stability in a slice of the wedge whose 
width appears to go to zero as /j —> -00. The problem here is that as /u —» -00 the DDE (1.2) is 
singularly perturbed and can be written as the so-called saw-tooth equation 

eu(t) = u{t) + Ku(t - a — cu(t)) 

where e — 1/fi and K = cr //r. This DDE had been studied in detail in [35] and for K > 1 
sufficiently large (corresponding to (j-i,cr) outside £*) the steady state is unstable, but there is 
an asymptotically stable slowly oscillating periodic solution. This periodic solution, known as 
the sawtooth solution, has unbounded gradient and a discontinuous profile in the singular limit. 
For parameter values inside the wedge the steady state is asymptotically stable, and for large 
and negative /./ there are no periodic solutions but a slowly decaying sawtooth-like oscillation 
can occur. Lyapunov-Razumikhin techniques based on bounding derivatives of solutions cannot 
perform well when those derivatives can be arbitrarily large. To improve the results in this case it 
would be necessary to define different sets £*^(<5, x) which take into account the structure of the 
oscillations and are hopefully much closer to E(k)(6, x) than the sets x) that we use here. 
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Figure 6: (a) The set X 2 \ Da which is part of the stability region of (1.2) for the constant delay case (c = 0) according to 
Bamea [1]. (b) The part of the stability domain outside Ea for the same problem found by Myshkis [37]. 


For p > 0 there is a significant improvement in the computed stability domain in going from 
k = 1 to k = 2 and a smaller improvement using k = 3. The largest value of p which satisfies 
P(1,0,1) ^ 1 can be computed from (5.1) which is quadratic in cr. Then non-negativity of the 
discriminant imposes the bound that p < (l/a)ln((l + V2)/2) w 0.1882 /a, as seen in Table 3. 

Although the parameter regions in which we can show asymptotic stability are independent of 
c, we will see in Section 6 that the basins of attraction do depend on c. 

6 . Basins of attraction 

Theorem 4.7 shows that for (p, cr) e U the ball 

[p : |M| < d 2 = ( y 1 r w+|,,1XB+|c|4 i>) (6.1) 

is contained in the basin of attraction of the steady-state of the state-dependent DDE (1.2) for 
k — 1,2 and 3. For fixed 6\ the radius of this ball gets smaller as k increases, but the value of 
ti] depends on k, p and cr, and some work is required to determine the largest such ball that is 
contained in the basin of attraction. In [31] we show that (6.1) can be improved when p < 0, so 
here we will consider (p, cr) e £ c , where cr < 0 ^ ji. Lemma 3.2 does not apply when /./ ^ 0, 
so there is no a priori bound on the solutions to (1.2) in this case. We present two examples 
which show that (1.2) can have unbounded solutions when p ^ 0, which also shows that the 
steady-state is not globally asymptotically stable when (p, cr) e Z r and gives an upper bound on 
the largest ball contained in its basin of attraction. For simplicity of exposition we suppose c > 0 
in this section, but the results can easily be extended to c < 0. We first consider p = 0. 

Example 6.1. Consider (1.2) with c>0, a > 0, p - 0 and cr e [-n/2a,-l/a) and for 6 e 
[-l/(ccr), a/c) let (fit) be Lipschitz continuous with 


<^>(0) = 6, ip(t) = -6 for t < -a - c6, 
ip{t ) € (-5, 6) for t e (-<7 - c6 ,0). 
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Figure 7: For fixed a = c = 1, /r = 0 and u £ \-nj2. 0], (a) supremum of <)| 6 (0, a/c] such that (p, cr) € { P(d \, c. k) < Oil 
with k = 1, 2, 3, and, (b) 62 = d\e *l‘ r K a+ l c l' y i) > 7 ’tie value of 6 * from Example 6.1 is also shown in both plots. 


Then while the deviated argument a(t, u(t)) = t - a — cu(t) ^ -a - c 6 we have 

u(t - a — cu(t )) = —6 and ii(t) — —crd. 

Hence, 

u(t) = h(l - crt) ^ 6 + for t ^ 0. (6.2) 

c 

But now ar(f, u(t)) ^ -a - c 6 for all t ^ 0 and (6.2) is valid for all t ^ 0. Thus for /j - 0, cr e 
[-n/ 2 a,-l/a) on the axis between the third and fourth quadrants of the stability region we have 
||<^|| = 6 and \u(t)\ —> 00 as t —> 00 . It follows that the steady state is not globally asymptotically 
stable and also that B( 0, 6 ) is not contained in its basin of attraction. Thus d* - -1 /ccr provides 
an upper bound on the radius of the largest ball contained in the basin of attraction. 

Figure 7 shows three bounds on the basin of attraction of the steady state of (1.2). The value 
of (V* from Example 6.1 gives an upper bound on the radius of the largest ball contained in the 
basin of attraction. Two lower bounds on the radius of the largest ball are also shown. The larger 
bound <5i gives the radius of the ball that Theorem 4.7 shows is contracted asymptotically to the 
steady state provided the solution is sufficiently differentiable. Lemma 2.3 is used to ensure that 
the solution remains bounded long enough to acquire sufficient regularity, and the growth in the 
solution allowed by that lemma results in the smaller radius 62 (as defined by (6.1)) of the ball 
that is contained in the basin of attraction for general continuous initial functions ip. We see that 
the bounds hi increase monotonically with k, but because of the exponential term in (6.1), the 
largest value of 62 is achieved with k = 1 in most of the interval for which F*(l, 0,1) < 1. 

Now consider the case of /r > 0. We can again derive an upper bound on the basin of attraction 
of the steady state when (jj. , cr) 6 

Example 6.2. Let a>0, c>0, /r>0 and (p., cr) e so cr < -/r < 0. Also let 

q{5) = -ap - 1 - ccr 6 - ln(c<5(yr - cr)). 

Note that q(l/(c(p - cr))) = -p(a + 1 /(p - cr)) < 0, while q'{5) = -ccr -1/6 < 0 for all 
6 e (0, l/(c(p - cr))). Also q{ 6 ) —> oo as 6 —* 0, hence there exists 6 * e (0,1 /(c(p - cr))) such 
that q( 6 *) = 0 and <5* is unique in this interval. 
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Figure 8: For fixed a = c = 1, cr = —1 and /./ > 0, (a) supremum of di € (0, a/c] such that (ji, cr) e { P(6\ , c, k) < <5i} with 
k = 2, 3, and, (b) 62 = Sie~ k ^ +,7 ^ a+ ^ Sl \ The value of <5* from Example 6.2 is also shown in (a). 


Suppose that the parameters are chosen so that 6 * < a/c. A sufficient (but not necessary) 
condition for this is cr < /r - 1 /a since this implies 1 /(c(ji - cr)) ^ a/c. Now let 6 e (<5*, a/c) so 
q( 6 ) < 0 and consider (1.2) with ip(t) Lipschitz continuous and 

<^(0) = 6 , ip(t ) = -6 for t < /iq{ 6 ), 

<p(t) e (- 6 , 6 ) for t € (jiq{ 6 ), 0), 


Then (1.2) has solution 


u(t ) = —- + 


/J - CT 


with u(t- a — cu{t)) — -6 for all t ^ 0. To see this note that 


(6.3) 


a(t, u(t)) = t - a - cu(t) — t - a - c6e^ - , 

H If*. 

with a(t, u(t)) —> -00 as t — > 00 . Differentiating the expression for a(t, u(t)) shows that 
a(t, u(t)) ^ nq(S) < 0 for all t ^ 0, with a(t , u(t)) = nq(S) when t — -[T l In (c6(fi - cr)). 
Hence, as in Example 6.1, we have ||^>|| = 6 and |w(f)| —> 00 as t —> 00 . The steady state is not 
globally asymptotically stable and the ball B{ 0, 6) is not contained in its basin of attraction. Thus 
6* provides an upper bound on the radius of the largest ball contained in the basin of attraction. 


For cr = — 1 and jj ^ 0, Figure 8 shows the same bounds 61,62 and 6 * on the radius of the basin 
of attraction of the steady state as were shown in Figure 7. Since these parameters are outside 
the set {P( 1,0, T) < 1) no bound is shown for k = 1. On nearly all of this interval k — 2 gives the 
largest lower bound 62 on the radius of a ball contained in the basin of attraction. 

Figure 9 shows these bounds on the basin of attraction in the cusp E c . The shaded region in 
Figure 9(c) denotes the portion of X, for which 6 * ^ a/c when a — c — 1, and hence 6 * from 
Example 6.2 gives an upper bound on the radius of the largest ball contained in the basin of 
attraction. The corresponding bounds 6 * are shown as contours within this region. Figure 9(a) 
and (b) shows the lower bounds 6 \ and 62 , along with the value of k that achieves the bound. 

In all three figures in this section, 62 is computed using (6.1) and t>i is obtained from solving 
P( 6 \, c, k ) = () 1 similarly to computations described in Section 5. 
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Figure 9: Plot of Z c for fixed a = c = 1 with contour plots of (a) the maximum 8\ e (0, a/c\ such that (ji,cr ) e 
{P(di, c, &) < di}, and (b) the 82 that maximizes 82 = (fe _ ^M+M)(0+kl<5 f or £ g (0, di]. Shading shows the value k e { 1,2, 3} 
for which the maximum is achieved, (c) The upper bound 8* from Example 6.2 for the radius of the largest ball B( 0, 8) 
contained in the basin of attraction of the zero solution to (1.2). 


7. Conclusions 

In this paper we have expanded upon the existing work on Lyapunov-Razumikhin techniques 
by providing results specifically tailored to DDEs with time-varying discrete delays including 
problems with state-dependent delays and vanishing delays. Our main results provide sufficient 
conditions for Lyapunov and asymptotic stability of steady state solutions of DDEs in Theo¬ 
rems 2.5 and 2.7 respectively. These conditions involve converting the DDE into a corresponding 
ODE problem with the delay terms treated as source terms that satisfy constraints. Our results 
require a Lipschitz condition on the right-hand side function / in ( 1. 1) instead of the more restric¬ 
tive Lipschitz condition on F in (1.4) required in Barnea [1], and do not require the construction 
of auxiliary functions as required by Hale and Verduyn Lunel [15]. Nevertheless we are able 
to show asymptotic stability, using a proof by contradiction showing that there cannot exist a 
solution which is not asymptotic to the steady state. 

We apply our results to the model state-dependent DDE ( 1 .2) in Sections 4-6. The main result 
of the application of Lyapunov-Razumikhin techniques to (1.2) is given as Theorems 4.7 where 
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we prove that the zero solution to (1.2) is asymptotically stable if (/u, <x) e {/ J ( 1,0, k) < 1), for 
k — 1, 2 or 3 and provide lower bounds on the basin of attraction. 

The parameter regions in which stability is proven in these theorems are compared in Section 5. 
As shown in Figure 5, the derived parameter regions grow as larger values of k are used, though 
for /r + 0, the derived stability region does not approach the entire known stability region X* as 
k —» oo (for reasons discussed in in Section 5). 

In Section 6 we consider (1.2) in the cusp X, where // > 0 and the steady state would be 
unstable without the delay term. In Examples 6.1 and 6.2 we constructed solutions which do not 
converge to the steady state for (ji, tr) e X, . These solutions provide us with an upper bound 6 * 
on the radius of the largest ball about the zero solution contained in the basin of attraction. In 
Figures 7-9 these upper bounds were compared with the lower bound 62 on the basin of attraction 
from (6.1). 

In the current work have studied stability through Lyapunov-Razumikhin techniques, but let 
us briefly compare and contrast this approach to the alternative, namely linearization. State- 
dependent DDEs have long been linearized by freezing the delays at their steady-state values and 
linearizing the resulting constant delay DDE [4, 5], This heuristic approach has recently been 
put on a rigorous footing. For a class of state-dependent DDEs which includes (1.2) with /j - 0, 
Gyori and Hartung [11] proved that the steady state of the state-dependent DDE is exponentially 
stable if and only if the steady state of the corresponding frozen-delay DDE is exponentially 
stable. In [12] they generalise this result to a class of nonautonomous problems which are linear 
except for the state-dependency. 

To compare and contrast our results with the linearization results of [12], we note that our 
results apply to a larger class of problems (1.1) than was considered in [12], and we prove both 
Lyapunov stability and asymptotic stability results, whereas [12] is concerned with exponen¬ 
tial stability. The results in [12] do apply directly to our model problem (1.2), and reveal the 
parameter region for which the steady state is exponentially stable. In contrast our Lyapunov- 
Razumikhin techniques are only able to deduce stability in part of this parameter region. 

Even though Lyapunov-Razumikhin techniques do not provide a proof of stability in the entire 
known stability region for (1.2), just as Lyapunov functions for ODEs do not always do so, they 
can nevertheless still be a very useful tool for studying stability in state-dependent DDEs. In 
particular our Lyapunov stability result is applicable to nonautonomous problems (for some of 
which rigorous linearization has yet to be derived) and the asymptotic stability result yield bounds 
on the basins of attraction which cannot be derived through linearization. 
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Appendix A. An explicit expression for the region P{6, c,2) < 6 


Here we prove Theorem 4.2. Let I\ and I 2 be defined by (4.12) and (4.13). Recall that J\ only 
applies when ^ r, in which case the integration does not have to be split into two intervals. 
For this case to occur we require u e [(rD\ - 1)<5, —pS/cr], which is only possible in the region 
where rD\ — 1 ^ -pier. When jyf < r the integration has to be broken into two parts and I2 
applies. In that case ii e | — 6 , min((/7L - I )d, -pd/crjj. We require the following lemmas. 

Lemma A.l. Let a > 0, c^O, cr ^ p and cr < -p. Let 6 e (0, \a/c\) and rD\ - 1 ^ -p/cr. If 
p ^ 0 then cr ^ — 1/r. If 11 < 0 then p e [(-3 + 2 V2)/r, 0] and 


cr > 


^(1 + /ur) + j -J 1 + 6 /ir + (fir) 2 




Proof. Let rD\ - 1 ^ -pier. Then 

r( sign(yu)yu - cr) - 1 = r(\p\ + |cr|) - 1 < r£>i - 1 < 

=> rcr + (1 - sign(/. i)pr)cr - p ^ 0. 

The boundary of the region where this inequality holds is 

cr = —7[^(1 - sign (p)pr) ± \ fl- sign(p)pr) 2 + 4pr]. 


(A.l) 


(A.2) 


(A.3) 


If p IP 0 then this simplifies to cr — 1 jr and cr — p. Since cr ' 0, and also p — cr — 0 satisfies 

(A.2) it follows that (A.l) holds for cr e [-1/r, 0] for the case p f 0. 

If At < 0 then (A. 3) simplifies to 

«■ = -^(i +pr)± j-y/l +6 pr + (pr) 2 ]. (A.4) 


Requiring 1 + 6 pr + (pr ) 2 ^ 0 yields pr e [-3 + 2 V2,0]. The lower bound on cr can be found by 
taking the lower boundary in (A.4) which attains its minimum at p — 0. This yields (A.l). □ 
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Lemma A.2. Let a > 0, cr ^ // and cr < —\jl. Let 6 e (0, |fl/c[). Define 

= fi£»i ln ( 1 -f^ r )- 1 ] ( 5. 

The following statements are true: 

(A) If rD\ - 1 > -///cr then the maximum of over u e [—c %—fj.d/cr] occurs at u„ if 

ft* < -[idler, and at —[id/cr otherwise. 

(B) If rD\ - 1 > -[i/cr and u * < —[id/cr then P( 6 , c, 2) ^ 5. 

(C) IfrDi - 1 ^ -/i/crthen sup fieHWrDj _ 1)tf] I 2 (u, 6 ) = I 2 ((rDi - 1)6,6) 

(D) IfrDx - 1 sj -fj/cr then sup- e[WJ| I\(u, 6) = J, (-[id/cr, 6) 

(E) IfrD\ ~ 1 ^ ~h/<t then l 2 ((rDi - l)b,(5) ^ I\(-/id/cr,d). 

Proof of (A). To find u, which maximises l 2 (u t , 6 ) consider the derivative 


6+u 


±I 2 (u,S) = er + f(S™- l). (A.5) 

At u — -6 this is positive. Setting the derivative equal to zero in (A.5) yields 

“* = [^1 ln(l-^ r )-l]h. 

Since 1 - ^e^ r e [0,1 ] if /i < 0 and 1 - ^e^ r > 1 if // > 0 then u„ > -6 in both cases. □ 

Proof of (B). First we need to prove that if ft, < -/id/cr and // > 0 then rl)\ — 1 ^ -// / cr. Let 
ii, < -[idler and /i > 0. Then jkl 2 (-[i 6 /cr, 6 ) < 0. Consider the term D\ ///, 


P P P \ P 


(A.6) 


Now consider the exponent of the second term in (A.5) with u = -/i 6 /cr, 

1 -P/o' = L _ P\P_ < 1 -P/^ = _P_ 

^ D i l cr) D\ ^ 1 -cr/n cr 


Thus, 


e , r + z {e -rhr _ 1) ^ fcM-fjS/o-,8) < 0. 

Isolating r in this expression yields r < ^ ln[^(l - Let x — / i/cr , then x e (—1,0) and 

(1 - e- x )/x > 1. Also, (1 - 1 /x) ln[(l - e~ x )/x] - 1 — x. These inequalities and (A.6) imply 

rDi- 1 < ^ in[2(1 - e~^)\ - 1 < (l - £)ln[z(l - e~^)\ - 1 ^ -£. 

Thus if rD\ - 1 > -///cr then ft, < -[idler can only occur if // < 0. 

6+u* 

Now let [i < 0 and ft* < -[idler. Then by setting JjJ 2 (ft,, d) - 0 in (A.5), ef* D v 5 - 1 = -^e / "'. 
Also, e Mr - cr/[i < 0 because // < 0 and cr /i < 0. Thus, 


> ~i) + ^ [-^ - - (t + ^ + i) 5 ^ r - 

But D\ > \/i\ + |cr[ > \ii/cr\\p\ + |cr[ = -(ji 2 /cr + cr) which implies y + ^ ^ < 

l 2 (fi*, 6 ) f 6 as required. 


0. Hence 

□ 
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Figure 10: Region where rD\ — 1 ^ -/i/cr is shown in brown and the region where e^ r + (o-/fi)(e^ r — 1) < 0 is shown in 
blue. These two regions do not intersect. 


Proof of (C). Let l 2 ({rD\ - 1)6,5) < 0. Then, e^ r + - 1) < 0, which can be rewritten as 


pe^"' 1 /-pr£^\ 

1 -e^ ~ r \ 1 - e^J ' 


(A.7) 


We show that the expression on the right-hand-side is continuous and decreases as p increases. 
Let x = iur, then 


lim -- f ZE?—\ = lim _i = -i, 

/j-> o r \ 1 — ev J x->o r \ \ — e x j r 


d 1 l -nre^ V 

c/p r \ 1 - e^ r / 


d xe x e x (l+x-e x ) 
dx 1 - e x (1 - e x ) 2 


So when p > 0, a necessary condition for (A. 7) to hold is cr < -1/r. From Lemma A.l, if 
rD\ - 1 ^ -p/cr and p > 0 then cr ^ -1/r. Thus Jrl 2 ((rDi - 1)5,6) ^ 0 if rD\ - 1 ^ -p/cr 
and p > 0. 

Now let rD\ - 1 ^ —y u/cr and p < 0. From Lemma A.l, pr e [-3 + 2 a/ 2,0] and cr ^ 
-i [i (1 + pr) - i Vl +6pr + (pr) 2 ]. Since =*£ ^ \ (1 + x) + \ Vl + 6* + x 2 for x 6 [-3 + 

2 V5,0], then cr ^ -4 )■ This contradicts (A. 7). Thus ^Ii((rD\ - 1)6,6) ^ 0 if rD\ - 1 ^ 

—p/cr and p < 0. Thus, 

r/L-l^-L => jj- u mrD x - 1)6,6) = 1) > 0. (A.8) 

This is shown in Figure 10. To finish the proof of (C), observe from (A.5) that 6) de¬ 
creases as u increases. Then by (A.8), sup Se[ _ tf (rfli _ 1)(5 ] 6) = l 2 ((rD\ - 1)5,5). □ 


Proof of (D). Let rD\ - 1 —p/cr. For all u G \(rD\ - 1)6, -pd/cr], 


|Ji(m,6) = - 1) = j-mrD, - 1)6,6). 

From (A.8), ^Ji(M) > 0 for all u e [-6, (rD x - 1)6]. Thus, sup fi€[(rW| _i^_ M/cr] J, (u, 6) = 
Ii(-/j.6/cr,6). □ 
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Proof of (E). From (4. 12) and (4.13) we find that 


h 6 ) - mrD x - 1)5,5) - (-£ - (rD x - 1)) 5 \e' lr + ^ - 1)]. 

For rD\ — 1 ^ —p/cr it follows from (A. 8) that this expression is non-negative. Thus, 
M-pd/a, 5) ^ J 2 ((rD, - 1)5,5). " □ 

Finally we can prove Theorem 4.2. 

Proof of Theorem 4.2. Note that this expression does not hold outside of {F(5, c, 2) < 5}. In order 
to prove this theorem, we need items (A)-(E) in Lemma A.2 which require Lemma A.l. 

Let rD\ - 1 > —p/cr. Then we can only have the two-part integration so I(u, 5, c, 2) = Ii(u, 5). 
From (A) and (B), P(6, c, 2) = J(-/r/cr, 5) if P(S, c, 2) < 5. 

Let rD\ - \ f —pier. Then we can have either the one-part or the two-part integration. 
From (C) and (D), I(u,6,c, 2) = max{J 2 ((rL>i - 1)5,5), Ii(-p6/cr,d)}. From (E), P(6,c, 2) = 
Ii(-p6/cr,6). □ 
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